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Abstract  (corn.) 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3-component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  96.7%  as  compared  to  99.6%  for  the  previous  reporting  period.  A 
total  of  2056  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  tested  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINESA  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  the  two  stations  in  Poland  have  been  recorded  and  pro¬ 
cessed  in  an  experimental  mode.  As  of  the  end  of  the  reporting  period  (30  Sep  1992)  data 
acquisition  from  these  two  stations  was  terminated,  in  accordance  with  the  terms  of  the 
contract.  Monthly  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis 
for  the  reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  involved  with  occasional  maintenance  of  equipment  for  FINESA 
and  work  in  connection  with  the  two  stations  in  Poland.  Other  activities  have  involved 
testing  of  the  NORSAR  communications  systems,  and  establishment  of  experimental 
small-aperture  arrays  at  sites  in  Spitsbergen  and  the  Kola  Peninsula  (see  Section  7.1). 

Starting  1  October  1991,  an  effort  has  begun  to  carry  out  a  complete  technical  refurbish¬ 
ment  of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  DARPA  and 
NTNF.  During  the  reporting  period,  efforts  have  focused  upon  evaluation  and  laboratory 
testing  of  technical  options  for  field  instrumentation,  in  particular  state-of-the-art  A/D 
converters,  data  acquisition  and  synchronization  devices.  During  the  next  few  months,  we 
plan  to  test  several  such  systems  under  realistic  operating  conditions  in  the  field.  Initial 
testing  of  some  systems  has  already  started.  When  these  studies  have  been  completed,  a 
recommendation  for  a  system  to  be  installed  will  be  presented  to  the  funding  agencies. 

Summaries  of  seven  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

Section  7.1  gives  a  technical  description  of  the  most  recent  extensions  of  the  Northern 
Europe  Regional  Array  Network.  Two  new  small-aperture  arrays  have  been  established; 
one  near  Apatity,  Russia,  and  one  on  the  Arctic  island  of  Spitsbergen.  The  Apatity  array 
was  installed  in  September  1992  as  part  of  an  agreement  on  scientific  cooperation  between 
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NORSAR  and  the  Kola  Science  Centre  of  the  Russian  Academy  of  Sciences.  This  9-ele¬ 
ment  array  comprises  a  center  site  and  two  concentric  rings,  and  has  an  aperture  of 
approximately  1  km.  All  sites  are  equipped  with  a  short-period  vertical  seismometer  of 
type  Geotech  S-500,  and  the  site  at  the  center  of  the  array  has  in  addition  two  horizontal 
seismometers  of  the  same  type.  Data  are  transmitted  to  Apatity  by  radio  link.  A  dedicated 
full  duplex  64  Kbps  satellite  link,  based  on  Norwegian  Telecom’s  NORSAT  B  system,  has 
been  installed  between  NORSAR  and  the  Kola  Regional  Seismology  Centre  in  Apatity.  A 
computer-to-computer  Ethernet  link  connects  the  two  data  centers. 

During  late  Octobet/early  November  1992  a  small- aperture  array  very  similar  to  the  Apat¬ 
ity  array  was  installed  on  the  island  of  Spitsbergen,  east  of  Longyearbyen.  Data  from  this 
array  are  transmitted  by  radio  and  terrestrial  link  to  Norwegian  Telecom’s  satellite  hub  sta¬ 
tion  at  Isfjord  Radio,  from  where  a  dedicated  simplex  64  Kbps  satellite  link  is  used  for 
transmission  of  the  continuous  data  to  Norway.  The  Spitsbergen  array  deployment  is  sup¬ 
ported  financially  by  Norwegian  oil  companies,  and  the  integration  of  this  array  into  the 
IMS  (Intelligent  Monitoring  System)  is  sponsored  by  DARPA. 

Section  7.2  presents  initial  processing  results  from  the  Apatity  array.  The  noise  level  is 
found  to  be  similar  to  that  of  ARCESS  at  frequencies  below  2  Hz,  and  slightly  higher  than 
ARCESS  at  higher  frequencies.  This  is  consistent  with  the  expectations.  Diurnal  noise 
variation  at  Apatity  is  similar  to  that  seen  at  NORESS.  The  Apatity  array  appears  to  pro¬ 
vide  excellent  noise  suppression  (Jn  or  better)  at  frequencies  above  2.5  Hz.  Initial  event 
processing  shows  that  frequency-wavenumber  processing  gives  well-defined  peaks  in  the 
F-K  diagram  both  for  P  and  S  type  phases.  A  particularly  noteworthy  feature  is  the  excel¬ 
lent  stability  in  azimuth  determination  of  local  events  using  the  low-frequency  Rg  phases. 
It  is  emphasized  that  these  results  are  preliminary  and  that  more  definite  conclusions  must 
await  analysis  of  more  extensive  data  sets. 

Section  7.3  discusses  the  use  of  regionalized  wave  propagation  characteristics  in  auto¬ 
matic  global  phase  association.  This  is  basically  an  attempt  to  establish  a  strategy  for 
retrieving  and  organizing  regionally  dependent  propagation  parameters  that  are  useful  for 
generalized  beamforming  of  global  network  data.  Examples  from  the  GSETT-2  data  base 
are  used  to  illustrate  features  such  as  the  distance-dependent  occurrence  of  seismic  phases, 
the  reliability  of  phase  labels  reported  by  NDCs  and  the  accuracy  of  one-station  event 
locations. 

Section  7.4  presents  and  compares  two  techiques  for  construction  of  a  uniform  grid  system 
covering  the  earth’s  surface.  The  first  technique  deploys  equidistant  grid  points  along 
equidistant  latitude  circles.  The  second  method  uses  triangulation  of  icosaeders.  The  sec¬ 
ond  technique  is  recommended  for  application  in  the  generalized  beamforming  process 
because  it  gives  more  effective  coverage  and  provides  a  well-defined  pointer  structure 
from  each  point  to  its  neighbors.  The  latter  property  is  important  for  beampacking,  since  it 
enables  us  to  refine  a  given  grid  in  a  straightforward  manner. 

Section  7.5  outlines  some  of  the  fundamental  concepts  for  using  the  generalized  beam- 
forming  (GBF)  method  to  conduct  phase  association  and  event  location  on  a  global  scale. 
Using  the  icosaeder-based  method  described  in  Section  7.4,  a  global  grid  of  beam-steering 
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points  is  defined.  Formulas  for  time  tolerances  and  slowness  vector  tolerances  are  devel¬ 
oped,  taking  into  account  both  the  tolerance  required  to  compensate  for  the  grid  spacing 
and  the  tolerances  due  to  effects  of  sampling  rates,  earth  inhomogeneities,  estimation 
uncertainties  and  other  random  errors.  A  step-by-step  description  for  a  global  GBF  algo¬ 
rithm  is  outlined.  An  example  of  processing  a  1-hour  time  interval  of  GSETT-2  detection 
data  is  presented.  This  interval  includes  a  small  earthquake  in  Tadjikistan,  and  this  earth¬ 
quake  is  correctly  processed  by  the  GBF  method  using  either  162, 642  or  2562  global  grid 
steering  points,  lhe  location  precision  is  illustrated  by  color  contour  maps,  and  is  natu¬ 
rally  highest  for  the  highest  point  density.  However,  it  is  noteworthy  that  an  acceptable 
solution  is  found  even  with  only  162  steering  points.  This  indicates  that  the  method  is  very 
robust  as  long  as  parameter  tolerances  are  adjusted  to  reflect  the  density  of  the  beam 
deployment. 

Section  7.6  reviews  available  NORSAR  detection  data  for  an  event  in  Ukraine  on  16  Sep 
1979,  which  was  recently  reported  in  Izvestiya  to  have  been  a  nuclear  explosion  of  1/3  kt 
yield,  detonated  next  to  a  Ukrainian  coal  mine.  Using  the  location  of  the  mining  town  as 
an  assumed  epicenter,  we  found  an  entry  in  the  NORSAR  automatic  detection  list  that  is 
consistent  with  an  origin  time  exactly  at  noon  (Moscow  time)  and  a  location  approxi¬ 
mately  as  given  in  the  press  report.  The  automatically  calculated  magnitude  was  mj,  =  3.3. 
Some  implications  of  this  case  study  for  seismic  monitoring  are  discussed. 

Section  7.7  is  a  review  of  induced  seismicity  caused  by  mining  activities  in  the  Khibiny 
Massif,  Kola  Peninsula.  It  is  demonstrated  that  the  extraction  of  large  volumes  of  rock 
mass  has  led  to  a  significant  increase  in  earthquakes  and  rockbursts  in  recent  years. 
Although  earthquakes  have  been  found  to  be  triggered  by  explosions  in  many  cases,  it  is 
too  early  to  give  a  definite  conclusion  abut  the  triggering  mechanism.  The  array  recently 
installed  near  Apatity  is  located  close  to  the  Khibiny  Massif,  and  will  provide  useful  data 
for  further  investigation  of  the  source  characteristics  of  earthquakes  and  mining  explo¬ 
sions  in  this  area. 
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1  Summary 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance  and  research 
activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Norwegian  Regional  Seismic 
Array  (NORESS)  and  the  Arctic  Regional  Seismic  Array  (ARCESS)  for  the  period 
1  April  -  30  September  1992.  Statistics  are  also  presented  for  additional  seismic  stations, 
which  through  cooperative  agreements  with  institutions  in  the  host  countries  provide  con¬ 
tinuous  data  to  the  NORSAR  Data  Processing  Center  (NPDC).  These  stations  comprise 
the  Finnish  Experimental  Seismic  Array  (FINESA),  the  German  Experimental  Seismic 
Array  (GERESS),  and  two  3-component  stations  in  Poland:  Ksiaz  and  Stary  Folwark. 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3-component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  96.7%  as  compared  to  99.6%  for  the  previous  reporting  period.  A 
total  of  2056  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  tested  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINESA  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  the  two  stations  in  Poland  have  been  recorded  and  pro¬ 
cessed  in  an  experimental  mode.  As  of  the  end  of  the  reporting  period  (30  Sep  1992)  data 
acquisition  from  these  two  stations  was  terminated,  in  accordance  with  the  terms  of  the 
contract.  Monthly  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis 
for  the  reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  involved  with  occasional  maintenance  of  equipment  fer  FINESA 
and  work  in  connection  with  the  two  stations  in  Poland.  Other  activities  have  involved 
testing  of  the  NORSAR  communications  systems,  and  establishment  of  experimental 
small -aperture  arrays  at  sites  in  Spitsbergen  and  the  Kola  Peninsula  (see  Section  7.1). 

Starting  1  October  1991,  an  effort  has  begun  to  cany  out  a  complete  technical  refurbish¬ 
ment  of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  DARPA  and 
NTNF.  During  the  reporting  period,  efforts  have  focused  upon  evaluation  and  laboratory 
testing  of  technical  options  for  field  instrumentation,  in  particular  state-of-the-art  A/D 
converters,  data  acquisition  and  synchronization  devices.  During  the  next  few  months,  we 
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plan  to  test  several  such  systems  under  realistic  operating  conditions  in  the  field.  Initial 
testing  of  some  systems  has  already  started.  When  these  studies  have  been  completed,  a 
recommendation  for  a  system  to  be  installed  will  be  presented  to  the  funding  agencies. 

Summaries  of  seven  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

Section  7.1  gives  a  technical  description  of  the  most  recent  extensions  of  the  Northern 
Europe  Regional  Array  Network.  Two  new  small-aperture  arrays  have  been  established; 
one  near  Apatity,  Russia,  and  one  on  the  Arctic  island  of  Spitsbergen.  The  Apatity  array 
was  installed  in  September  1992  as  part  of  an  agreement  on  scientific  cooperation  between 
NORSAR  and  the  Kola  Science  Centre  of  the  Russian  Academy  of  Sciences.  This  9-ele¬ 
ment  array  comprises  a  center  site  and  two  concentric  rings,  and  has  an  aperture  of 
approximately  1  km.  All  sites  are  equipped  with  a  short-period  vertical  seismometer  of 
type  Geotech  S-500,  and  the  site  at  the  center  of  the  array  has  in  addition  two  horizontal 
seismometers  of  the  same  type.  Data  are  transmitted  to  Apatity  by  radio  link.  A  dedicated 
full  duplex  64  Kbps  satellite  link,  based  on  Norwegian  Telecom’s  NORSAT  B  system,  has 
been  installed  between  NORSAR  and  the  Kola  Regional  Seismology  Centre  in  Apatity.  A 
computer-to-computer  Ethernet  link  connects  the  two  data  centers. 

During  late  Octobei/early  November  1992  a  small- aperture  array  very  similar  to  the  Apat¬ 
ity  array  was  installed  on  the  island  of  Spitsbergen,  east  of  Longyearbyen.  Data  from  this 
array  are  transmitted  by  radio  and  terrestrial  link  to  Norwegian  Telecom's  satellite  hub  sta¬ 
tion  at  Isfjord  Radio,  from  where  a  dedicated  simplex  64  Kbps  satellite  link  is  used  for 
transmission  of  the  continuous  data  to  Norway.  The  Spitsbergen  array  deployment  is  sup¬ 
ported  financially  by  Norwegian  oil  companies,  and  the  integration  of  this  array  into  the 
IMS  (Intelligent  Monitoring  System)  is  sponsored  by  DARPA. 

Section  7.2  presents  initial  processing  results  from  the  Apatity  array.  The  noise  level  is 
found  to  be  similar  to  that  of  ARCESS  at  frequencies  below  2  Hz,  and  slightly  higher  than 
ARCESS  at  higher  frequencies.  This  is  consistent  with  the  expectations.  Diurnal  noise 
variation  at  Apatity  is  similar  to  that  seen  at  NORESS.  The  Apatity  array  appears  to  pro¬ 
vide  excellent  noise  suppression  (Jn  or  better)  at  frequencies  above  2.5  Hz.  Initial  event 
processing  shows  that  frequency-wavenumber  processing  gives  well-defined  peaks  in  the 
F-K  diagram  both  for  P  and  S  type  phases.  A  particularly  noteworthy  feature  is  the  excel¬ 
lent  stability  in  azimuth  determination  of  local  events  using  the  low-frequency  Rg  phases. 
It  is  emphasized  that  these  results  are  preliminary  and  that  more  definite  conclusions  must 
await  analysis  of  more  extensive  data  sets. 

Section  7.3  discusses  the  use  of  regionalized  wave  propagation  characteristics  in  auto¬ 
matic  global  phase  association.  This  is  basically  an  attempt  to  establish  a  strategy  for 
retrieving  and  organizing  regionally  dependent  propagation  parameters  that  are  useful  for 
generalized  beamforming  of  global  network  data.  Examples  from  the  GSETT-2  data  base 
are  used  to  illustrate  features  such  as  the  distance-dependent  occurrence  of  seismic  phases, 
the  reliability  of  phase  labels  reported  by  NDCs  and  the  accuracy  of  one-station  event 
locations. 
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Section  7.4  presents  and  compares  two  techiques  for  construction  of  a  uniform  grid  system 
covering  the  earth's  surface.  The  first  technique  deploys  equidistant  grid  points  along 
equidistant  latitude  circles.  The  second  method  uses  triangulation  of  icosaeders.  The  sec¬ 
ond  technique  is  recommended  for  application  in  the  generalized  beamforming  process 
because  it  gives  more  effective  coverage  and  provides  a  well-defined  pointer  structure 
from  each  tXjint  to  its  neighbors.  The  latter  property  is  important  for  beampacking,  since  it 
enables  us  to  refine  a  given  grid  in  a  straightforward  manner. 

Section  7.5  outlines  some  of  the  fundamental  concepts  for  using  the  generalized  beam¬ 
forming  (GBF)  method  to  conduct  phase  association  and  event  location  on  a  global  scale. 
Using  the  icosaeder-based  method  described  in  Section  7.4,  a  global  grid  of  beam-steering 
points  is  defined.  Formulas  for  time  tolerances  and  slowness  vector  tolerances  are  devel¬ 
oped,  taking  into  account  both  the  tolerance  required  to  compensate  for  the  grid  spacing 
and  the  tolerances  due  to  effects  of  sampling  rates,  earth  inhomogeneities,  estimation 
uncertainties  and  other  random  errors.  A  step-by-step  description  for  a  global  GBF  algo¬ 
rithm  is  outlined.  An  example  of  processing  a  1-hour  time  interval  of  GSETT-2  detection 
data  is  presented.  This  interval  includes  a  small  earthquake  in  Tadjikistan,  and  this  earth¬ 
quake  is  correctly  processed  by  the  GBF  method  using  either  162, 642  or  2562  global  grid 
steering  points.  The  precision  is  illustrated  by  color  contour  mapos,  and  is  naturally  high¬ 
est  for  the  highest  point  density.  However,  it  is  noteworthy  that  an  acceptable  solution  is 
found  even  with  only  162  steering  points.  This  indicates  that  the  method  is  very  robust  as 
long  as  parameter  tolerances  are  adjusted  to  reflect  the  density  of  the  beam  deployment. 

Section  7.6  reviews  available  NORSAR  detection  data  for  an  event  in  Ukraine  on  16  Sep 
1979,  which  was  recently  reported  in  Izvestiya  to  have  been  a  nuclear  explosion  of  1/3  kt 
yield,  detonated  next  to  a  Ukrainian  coal  mine.  Using  the  location  of  the  mining  town  as 
an  assumed  epicenter,  we  found  an  entry  in  the  NORSAR  automatic  detection  list  that  is 
consistent  with  an  origin  time  exactly  at  noon  (Moscow  time)  and  a  location  approxi¬ 
mately  as  given  in  the  press  report.  The  automatically  calculated  magnitude  was  m(,  =  3.3. 
Some  implications  of  this  case  study  for  seismic  monitoring  are  discussed. 

Section  7.7  is  a  review  of  induced  seismicity  caused  by  mining  activities  in  the  Khibiny 
Massif,  Kola  Peninsula.  It  is  demonstrated  that  the  extraction  of  large  volumes  of  rock 
mass  has  led  to  a  significant  increase  in  earthquakes  and  rockbursts  in  recent  years. 
Although  earthquakes  have  been  found  to  be  triggered  by  explosions  in  many  cases,  it  is 
too  early  to  give  a  definite  conclusion  abut  the  triggering  mechanism.  The  array  recently 
installed  near  Apatity  is  located  close  to  the  Khibiny  Massif,  and  will  provide  useful  data 
for  further  investigation  of  the  source  characteristics  of  earthquakes  and  mining  explo¬ 
sions  in  this  area. 
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2  NORSAR  Operation 

2.1  Detection  Processor  (DP)  operation 

There  have  been  88  breaks  in  the  otherwise  continuous  operation  of  the  NORSAR  online 
system  within  the  6-month  reporting  interval.  The  uptime  percentage  for  the  period  is 
96.1%  as  compared  to  99.6%  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  downtime  for  the 
days  between  1  April  and  30  September  1992.  The  monthly  recording  times  and  percent¬ 
ages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  43 

b)  Stops  related  to  program  work  or  error  0 

c)  Hardware  maintenance  stops  8 

d)  Power  jumps  and  breaks  3 

e)  TOD  error  correction  0 

0  Communication  lines  34 


The  total  downtime  for  the  period  was  142  hours  and  14  minutes.  The  mean-time  - 
between-failures  (MTBF)  was  2.0  days,  as  compared  to  4.6  for  the  previous  period. 


J.  Torstveit 
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Fig.  2.1.1.  Detection  Processor  uptime  for  July  (top),  August  (middle)  and  September 
(bottom)  1992. 
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Date  Time  Cause 


04  Apr 

0057 

-  1512 

Hardware  failure 

06  Apr 

0930 

-  0956 

Hardware  maintenance 

13  Apr 

2333 

- 

Hardware  failure 

14  Apr 

-  1511 

Hardware  failure 

21  Apr 

2052 

-  2140 

Hardware  failure 

23  Apr 

0653 

-  0841 

Hardware  maintenance 

27  Apr 

0540 

-  0648 

Line  failure 

28  Apr 

1826 

-  1907 

Hardware  failure 

14  May 

0905 

-  1034 

Line  failure 

24  May 

1125 

-  1358 

Aircondition  failure 

01  Jun 

2200 

- 

Hardware  failure  after  power  break 

02  Jun 

-  1352 

Hardware  failure  after  power  break 

05  Jun 

0223 

-  0854 

Hardware  failure 

07  Jun 

0003 

- 1319 

Power  break 

08  Jun 

1153 

-  1445 

Hardware  failure  due  to  thunderstorm 

10  Jun 

2059 

-  2158 

Hardware  failure 

22  Jun 

0712 

-  0810 

Hardware  maintenance 

23  Jun 

0200 

-  0311 

Hardware  failure 

27  Jun 

0636 

-  0711 

Hardware  maintenance 

14  Jul 

1123 

-  1158 

Hardware  failure 

18  Jul 

1800 

- 

Hardware  failure 

19  Jul 

-  1128 

Hardware  failure 

29  Jul 

1423 

-  1514 

Hardware  failure 

01  Aug 

0120 

-  0200 

Hardware  failure 

07  Aug 

1914 

-  2000 

Hardware  failure 

17  Aug 

1602 

-  1741 

Hardware  failure 

28  Aug 

2352 

- 

Aircondition  failure 

29  Aug 

-  0816 

Aircondition  failure 

29  Aug 

2024 

-2107 

Hardware  failure 

30  Aug 

1643 

-  1728 

Hardware  failure 

31  Aug 

0257 

-  0341 

Hardware  failure 

31  Aug 

1434 

-  1738 

Hardware  failure 

31  Aug 

2212 

-  2332 

Hardware  failure 

01  Sep 

0005 

-  0440 

Hardware  failure 

01  Sep 

0715 

-  0755 

Hardware  maintenance 

01  Sep 

1136 

-  1240 

Hardware  maintenance 

02  Sep 

0233 

-  1215 

Hardware  failure 

19  Sep 

0119 

-  0203 

Hardware  failure 

22  Sep 

0359 

-  0451 

Hardware  failure 

30  Sep 

1747 

-  1841 

Hardware  failure 

Table  2.1.1.  The  major  downtimes  in  the  period  1  April  -  30  September  1992. 


7 


NORSAR  Sci.  Rep.  1-92/93 


November  1992 


Month  DP  Uptime 

DP  Uptime 

No.  of 

No.  of  Days 

DP  MTBF* 

Hours 

% 

DP  Breaks 

with  Breaks 

(days) 

Apr  92 

684.38 

95.09 

16 

11 

1.7 

May  92 

739.33 

99.40 

6 

6 

4.4 

Jun  92 

676.02 

93.89 

22 

16 

1.2 

Jul  92 

723.57 

97.31 

13 

10 

2.2 

Aug  92 

724.51 

97.31 

20 

13 

1.4 

Sep  92 

700.45 

97.34 

11 

7 

2.4 

96.74 

88 

63 

2.0 

♦Mean-time-between-failures  =  total  uptime/no.  of  up  intervals. 


Table  2.1.2.  Online  system  performance,  1  April  -  30  September  1992. 
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2.2  Array  communications 

General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  throughout  the  report¬ 
ing  period.  Also  this  period  single,  groups  and  all  systems  (simultaneously)  have  been 
affected.  Types  of  events  have  varied,  as: 

•  Bad  communications  cable 

•  Reduced  line  level 

•  SLEM  (stuck) 

•  SLEM  (power) 

•  Power  (SA) 

•  Power  (NDPC) 

•  Sync  problems 

•  Modems  (CTV) 

•  Modem  (Avanti  230  Kbit  NDPC)  caused  bad  communications  between  Morlcomp 
and  2701  communications  adapter 

•  NTA  carrier 

•  Lightning 

•  IBM  disk  storage 

•  Defective  UPS  batteries  (NDPC) 

Detailed  Summary 

April  (weeks  14-18),  303-3.5.92 

Reliable  performance  for  all  systems  weeks  1 4  and  15,  except  for  06C,  which  was  affected 
3  and  4  April.  Besides  reliable  operation  01 B,  02B,  02C  and  06C  week  17.  All  systems 
were  affected  weeks  16  and  18.  01A,  03C  and  04C  were  affected  week  17.  An  IBM  disc 
storage  device,  the  NTA  carrier  system  and  the  Modcomp  have  caused  most  of  the  prob¬ 
lems.  Average  outage  in  April,  individual  weeks  (5): 

Week  14  (-06C)  :  0.0005%  0.115%  (all) 

Week  15  (all)  :  0.0006% 

Week  16  (all)  :  9.821% 

Week  17  (-01A.03C.04C)  :  0.001%  16.96%  (all) 

Week  18  (all)  ;  15.17% 

May  (weeks  19-22),  4-31.5.92 

The  NTA  carrier  frequency  system  failed  again  between  3  and  4  May,  week  21.  All  sys¬ 
tems  were  affected.  A  line  test  was  carried  out  between  NDPC  and  03C  14  May  in  cooper¬ 
ation  with  NMC/Hamar.  No  errors  were  observed  in  C-loop  (digital  loop).  02B  was  down 
between  20  and  26  May  due  to  a  power  outage.  02C  was  affected  between  22  and  24  May, 
probably  caused  by  the  line.  01B  was  down  between  26  May  and  1  June.  Average  outage 
in  May,  individual  weeks: 
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Week  19  (all)  :  4.169% 

Week  20  (all)  :  0.0008% 

Week  21  (-02C)  :  0.002% 

Week  22  (-01B,02B)  :  0.0005% 


June  (weeks  23-26),  1-28.6.92 

Without  warning  the  local  power  company  disconnected  the  power  1  June  at  2200  hrs. 
The  UPS  (Uninterrupted  Power  Supply)  failed  and  therefore  all  systems  lost  power  imme¬ 
diately.  The  alarm  system  did  not  warn  the  person  on  duty,  and  as  a  result  of  this  the  sys¬ 
tems  were  not  restarted  until  the  next  day  at  1355  hrs. 

01 B  went  down  4  June  due  to  a  faulty  SLEM  power  unit  (which  also  supplies  the  modem), 
caused  excessive  amounts  of  communications  errors.  12  June  the  power  unit  was  replaced. 

5  June  NORSAR  was  down  approximately  6  hrs  30  min  in  connection  with  replacing  a 
power  supply  in  a  disc  controller. 

7  June  a  short  power  outage  caused  vital  equipment,  such  as  the  2701  adapter  and  a  disc 
controller  to  drop  power.  In  this  connection  the  communications  systems  was  down  for 
about  13  hrs. 

9  June  02C  went  down  20.38  hrs.  A  Modcomp  restart  the  next  day  restored  the  subarray 
operation. 

11  June  a  01B  communications  test  revealed  bad  C-loop  (digital  loop).  12  June  01B 
SLEM  power  unit  was  replaced  and  normal  operation  was  restored. 

16  June  A  comtest  on  the  01 A  B/C-Ioop  revealed  a  bad  communications  cable.  NTA/ 
Hamar  tried  to  improve  the  conditions  by  using  two  other  pairs  in  the  cable  but  without 
success.  A  new  communications  test  22  June  only  confirmed  the  degraded  quality. 
According  to  NTA/Hamar  the  cable  repair  had  to  be  postponed  because  the  damage  was 
located  to  that  part  of  the  cable  crossing  a  com  field. 

24  June  at  2226  hrs  spikes  and  “time  mismatch”  was  observed  on  the  NORSAR  system. 

25  June  the  Modcomp  was  restarted  and  normal  operation  restored. 

Average  outage  in  June,  individual  weeks: 


Week  23  (-01B)  :  0.003% 

Week  24  (-01B.02C)  :  0.021% 

Week  25  (all)  :  0.035% 

Week  26  (-01A)  :  0.004% 


July  (weeks  27-31),  29.6-2.8.92 

Also  in  July  the  NORSAR  communications  systems  were  affected  either  individually  or 
all  simultaneously.  The  individually  affected  systems  were  caused  by  bad  communications 
cables,  short  outages,  lightning  in  the  subarray  area,  sync  problems,  etc.  Simultaneously 
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affected  systems  were  probably  caused  by  shot  interruptions  in  the  NTA  group  transmis¬ 
sion  equipment.  Also  the  NDPC  communications-related  equipment  has  been  considered 
as  a  possible  source. 

01A  was  inoperative  also  throughout  July.  02C  and  04C  were  down  between  1  and  2  July. 
02C  was  affected  again  6  and  7  July.  All  systems  were  affected  between  18  and  19  July. 
The  Modcomp  was  restarted  and  the  systems  resumed  operation  (-01A,01B,02B).  04C 
and  06C  were  affected  14  July.  19  July  02 C,  03C,  04C  and  06C  were  affected  for  approxi¬ 
mately  1  hour.  20  July  01B  resumed  operation;  02B  on  21  July  (both  had  been  down  since 
19  July).  02B  was  affected  again  27-28  July  caused  by  lightning.  02C  indicated  sync  prob¬ 
lems  29  July  between  0000  and  0300  hrs. 

In  the  period  the  Modcomp  was  restarted  16  times  in  order  to  restore  the  operation  of  one 
or  several  subarrays,  but  also  in  order  to  align  TOD  band  NORSAR  data. 

Average  outages  in  July,  individual  weeks: 


Week  27  (-01A)  0.007% 

Week  28  (-01A.02C)  0.009% 

Week  29  (N/A) 

Week  30  (-01A,01B,02B)  0.009% 

Week  31  (-01  A)  0.043% 


August  (weeks  32-35),  3-30.8.92 

01A  was  also  out  of  operation  also  in  August.  01B,  02B,  04C  and  06C  have  been  affected 
this  period. 

01B  performance  was  reduced  9-10  August,  but  the  SLEM  unit  was  replaced  and  opera¬ 
tion  established  again  20  August. 

02B  was  affected  5-7  and  11-13  August  in  connection  with  lightning  in  the  02B  area.  04C 
was  affected  by  a  thunderstorm  12-13  August,  and  06C  was  down  between  10  and  11 
August,  probably  due  to  short  line  outages. 

All  system  (-02C)  dropped  out  10  August  between  20  and  23  hrs.  After  a  Modcomp  restart 
0453  hrs  1 1  August  the  systems  resumed  operation. 

Average  outages  in  August,  individual  weeks  (4): 

Week  32  (-01A.02B)  0.016% 

Week  33  (N/A) 

Week  34  (-01A.01B)  0.002% 

Week  35  (-01A)  0.009% 
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September  (weeks  36-39),  31.8-27.9.92 

1  September  between  1136  and  1240  hrs  the  communication  between  the  IBM  2701  and 
the  Modcomp  was  broken,  caused  by  a  faulty  Avanti  230  Kbit  modem.  After  repair,  the 
data  exchange  started  again. 

2  September  a  disc  failed  which  again  caused  loss  of  NORSAR  data  from  0233  hrs  until 
the  disc  was  replaced  1215  hrs. 

06C  was  affected  2  September  between  10  and  13  hrs.  A  Modcomp  restart  reinitiated  the 
system.  01B  was  affected  22-23  September  caused  by  a  broken  communications  cable 
between  Kjeller  and  Lillestrom.  Power  at  02B  was  broken  a  few  hours  23  September,  but 
did  not  affect  the  subarray  operation  due  to  satisfactory  operation  of  the  backup  batteries. 

Average  outages  in  September,  individual  weeks  (4): 


Week  36  (-01  A)  :  0.028% 

Week  37  (-01A)  :  0.0008% 

Week  38  (-01  A)  :  0.012% 

Week  39  (-01  A)  :  0.298% 

(-01A,  01B)  0.0007% 


O.A.  Hansen 
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2.3  NORSAR  Event  Detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  lists  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detec¬ 
tor,  the  total  .lumber  of  detections  processed  by  the  automatic  event  processor  (EPX)  and 
the  total  number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases 
and  total). 


Total 

DPX 

Total 

EPX 

P-phases 

Accepted  Phases 
Core  Phases 

Sum 

Daily 

Apr  92 

10400 

1308 

245 

64 

309 

10.3 

May  92 

6100 

1301 

384 

104 

488 

15.7 

Jun  92 

7600 

1301 

296 

87 

383 

12.8 

Jul  92 

9825 

1505 

359 

69 

428 

13.8 

Aug  92 

9225 

1402 

383 

62 

445 

14.4 

Sep  92 

9000 

1188 

269 

53 

322 

10.7 

1936 

439 

2375 

13.0 

Table  23.1.  Detection  and  Event  Processor  statistics,  1  April  -  30  September  1992. 


NORSAR  Detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  092, 
through  day  274, 1992,  was  52,129,  giving  an  average  of  285  detections  per  processed  day 
(183  days  processed).  Table  2.3.2  shows  daily  and  hourly  distribution  of  detections  for 
NORSAR. 

B.  Paulsen  Gammelby 
T.  Scheyen 


.  DPX  Hourly  distribution  of  detections 
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Table  2.3.2.  Daily  and  hourly  distribution  of  NORSAR  detections.  For  each  day  is  shown  number  of  detections  within  each 
hour  of  the  day  and  number  of  detections  for  that  day.  The  end  statistics  give  total  number  of  detections  distributed  for  each 
hour  and  the  total  sum  of  detections  during  the  period.  The  averages  show  number  of  processed  days,  hourly  distribution  and 
average  per  processed  day. 
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3  Operation  of  regional  arrays 

3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 

Table  3.1.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  98.55%  as  compared  to  99.78%  during  the  previous 
reporting  period. 


Date  Time  Cause 


09  Apr 

0127 

-  0225 

Transmission  line  failure 

09  Apr 

1007 

-  1426 

Transmission  line  failure 

05  May 

2342 

- 

Transmission  line  failure 

06  May 

-  0008 

Transmission  line  failure 

01  Jun 

2200 

- 

Power  break  at  NDPC 

02  Jun 

-  0454 

Power  break  at  NDPC 

02  Jun 

0454 

- 

Reduced  data  quality  (HUB  problems) 

03  Jun 

-  0130 

Reduced  data  quality  (HUB  problems) 

03  Jun 

1207 

- 

Reduced  data  quality  (HUB  problems) 

04  Jun 

-  0100 

Reduced  data  quality  (HUB  problems) 

04  Jun 

1130 

-  2200 

Reduced  data  quality  (HUB  problems) 

11  Aug 

1607 

-  1723 

Transmission  line  maintenance 

11  Aug 

1731 

-  1923 

Transmission  line  maintenance 

11  Aug 

1927 

-  2323 

Transmission  line  maintenance 

11  Aug 

2329 

- 

Transmission  line  maintenance 

12  Aug 

-  0123 

Transmission  line  maintenance 

12  Aug 

0145 

-  0324 

Transmission  line  maintenance 

12  Aug 

1108 

-  1124 

Transmission  line  maintenance 

12  Aug 

1151 

-  1324 

Transmission  line  maintenance 

27  Sep 

0059 

-  0200 

Transmission  line  failure 

Table  3.1.1.  Interruptions  in  recording  of  NORESS  data  at  NDPC,  1  April  -  30  September 
1992. 


Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 
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April 

99.20 

May 

99.88 

June 

96.23 

July 

99.95 

August 

96.17 

September 

99.86 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  3.1.1.  NORESS  data  recording  uptime  for  April  (top),  May  (middle)  and  June  (bot¬ 
tom)  1992. 
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Fig.  3.1.1.  (corn.)  NORESS  data  recording  uptime  for  July  (top),  August  (middle)  and 
September  (bottom)  1992. 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjeller 

Table  3.2.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  99.25  %  as  compared  to  97.28%  for  the  previous  report¬ 
ing  period.. 


Date  Time  Cause 


01 

Apr 

1104 

-  1130 

26 

Apr 

1718 

-  2005 

28 

Apr 

1245 

-  2106 

11 

May 

0926 

-  0936 

01 

Jun 

2200 

- 

02 

Jun 

-  0457 

06  Jun 

2040 

-  2130 

21 

Jun 

1826 

-  2031 

22 

Jun 

1744 

-  1754 

14 

Jul 

2049 

-  2137 

17 

Aug 

1658 

-  1710 

17 

Aug 

1733 

-  1748 

17 

Aug 

1902 

-  1917 

07 

Sep 

0646 

-  0656 

07 

Sep 

2310 

- 

08 

Sep 

-  0052 

Service  on  the  satellite  modem 
Power  break  HUB 
Hardware  failure  NDPC 
Service  on  the  satellite  modem 
Power  break  NDPC 
Power  break  NDPC 
Hardware  failure  NDPC 
Service  HUB 
Service  HUB 
Unknown 

Satellite  link  failure 
Satellite  link  failure 
Satellite  link  failure 
Hardware  service  NDPC 
Satellite  link  failure 
Satellite  link  failure 


Table  3.2.1.  The  main  interruptions  in  recording  of  ARCESS  data  at  NDPC,  1  April  -  30 
September  1992. 

Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

97.74% 

May 

99.97% 

June 

98.37% 

July 

99.86% 

August 

99.84% 

September 

99.71% 

Fig.  3.2.1 .  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  April  (top).  May  (middle)  and  June(bot- 
tom)  1992. 
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3.3  Recording  of  FINESA  data  at  NDPC,  Kjelier 

The  average  recording  time  was  93.75%  as  compared  to  95.5%  for  the  previous  period.  As 
can  be  seen  from  Table  3.3.1  below,  the  main  reason  for  the  downtime  is  transmission  line 
failure  and  Hub  failure. 


Date  Time  Cause 


23  Apr 

2328 

- 

Transmission  line  failure 

24  Apr 

-  0518 

Transmission  line  failure 

09  May 

1042 

- 

Transmission  line  failure 

11  May 

-  0917 

Transmission  line  failure 

22  May 

1000 

-  1021 

Transmission  line  failure 

01  Jun 

2200 

- 

Power  failure  NDPC 

02  Jun 

-  0458 

Power  failure  NDPC 

18  Jun 

1833 

-  2140 

Transmission  line  failure 

20  Jun 

0853 

- 

Transmission  line  failure 

22  Jun 

-  1310 

Transmission  line  failure 

05  Jul 

0309 

- 

Transmission  line  failure 

06  Jul 

-  0727 

Transmission  line  failure 

13  Jul 

1359 

- 

Transmission  line  failure 

14  Jul 

-  0815 

Transmission  line  failure 

27  Jul 

1116 

-  1143 

Transmission  line  failure 

31  Jul 

1240 

- 

Hardware  failure  HUB 

03  Aug 

-  1218 

Hardware  failure  HUB 

06  Aug 

1137 

-  1149 

Transmission  line  failure 

10  Aug 

2038 

-  2109 

Transmission  line  failure 

10  Aug 

2119 

- 

Transmission  line  failure 

11  Aug 

-  0808 

Transmission  line  failure 

17  Aug 

1229 

-  1327 

Transmission  line  failure 

01  Sep 

0335 

- 

Transmission  line  failure 

02  Sep 

-  0638 

Transmission  line  failure 

Table  33.1.  The  main  interruptions  in  recording  of  FINESA  data  at  NDPC,  1  April  -  30 
September  1992. 


Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmission  lines,  data  center  operation)  affecting  this  task  were 
as  follows: 
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April 

99.13% 

May 

93.65% 

June 

91.25% 

July 

92.09% 

August 

90.17% 

September 

96.22% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
FINESA  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J,  Torstveit 
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Fig.  3.3.1.  FINESA  data  recording  uptime  for  April  (top),  May  (middle)  and  June  (bot¬ 
tom)  1992. 
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Fig.  3.3.1.  FINESA  data  recording  uptime  for  July  (top),  August  (middle)  and  September 
(bottom)  1992. 
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3.4  Event  detection  operation 

This  section  reports  results  from  one-array  automatic  processing  using  signal  processing 
recipes  and  “ronapp”  recipes  for  the  ep  program  (NORSAR  Sci.  Rep.  No  2-8&89). 

Three  systems  are  in  parallel  operation  to  associate  detected  phases  and  locate  events: 

1.  The  ep  program  with  “ronapp”  recipes  is  operated  independently  on  each  array  to 
obtain  simple  one-array  automatic  solutions. 

2.  The  Generalized  Beamforming  method  (GBF)  (see  F.  Ringdal  and  T.  Kvaema 
(1989),  A  mulitchannel  processing  approach  to  real  time  network  detection,  phase 
association  and  threshold  monitoring,  BSSA  Vol  79,  no  6, 1927-1940)  processes  the 
four  arrays  jointly  and  presents  locations  of  regional  events. 

3.  The  IMS  system  is  operated  on  the  same  set  of  arrivals  as  ep  and  GBF  and  reports 
also  teleseismic  events  in  addition  to  regional  ones. 

IMS  results  are  reported  in  section  3.5  and  GBF  results  in  section  3.6. 

In  addition  to  these  three  event  association  processes,  we  are  running  test  versions  of  the 
so-called  Threshold  Monitoring  (TM)  process.  This  is  a  process  that  monitors  the  seismic 
amplitude  level  at  the  four  regional  arrays  continuously  in  time  to  estimate  the  upper  mag¬ 
nitude  limit  of  an  event  that  might  go  undetected  by  the  network.  The  current  TM  process 
is  beamed  to  several  sites  of  interest,  including  the  Novaya  Zemlya  test  site.  Simple  dis¬ 
plays  of  so-called  threshold  curves  reveal  instants  of  particular  interest;  i.e.,  instants  when 
events  above  a  certain  magnitude  threshold  may  have  occurred  in  the  target  region. 
Results  from  the  three  processes  described  above  are  used  to  help  resolve  what  actually 
happened  during  these  instances. 

NORESS  detections 

The  number  of  detections  (phases)  reported  from  day  092, 1992,  through  day  274, 1992, 
was  42,587,  giving  an  average  of  233  detections  per  processed  day  (183  days  processed). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  092, 1992,  through  274, 1992, 2556  local  and  regional  events  were  located  by 
NORESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  14.0  events  per  processed  day  (183  days  processed).  60%  of  these  events  are  within 
300  km,  and  85%  of  these  events  are  within  1000  km. 

ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  092, 1992,  through  day  274, 1992, 
was  92,044,  giving  an  average  of  503  detections  per  processed  day  (182  days  processed). 
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ARC  ESS  detections 

The  number  of  detections  (phases)  reported  during  day  092, 1992,  through  day  274,  1992, 
was  92,044,  giving  an  average  of  503  detections  per  processed  day  (182  days  processed). 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  092, 1992,  through  274, 1992,  3875  local  and  regional  events  were  located  by 
ARCESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
21.3  events  per  processed  day  (182  days  processed).  47%  of  these  events  are  within 
300  km,  and  86%  of  these  events  are  within  1000  km. 

FINESA  detections 

The  number  of  detections  (phases)  reported  during  day  092, 1992,  through  day  274, 1992, 
was  43,036,  giving  an  average  of  240  detections  per  processed  day  (179  days  processed). 

Table  3.4.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 

Events  automatically  located  by  FINESA 

During  days  092, 1992,  through  274, 1992,  2973  local  and  regional  events  were  located  by 
FINESA,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  16.6  events  per  processed  day  (179  days  processed).  71%  of  these  events  are  within 
300  km,  and  88%  of  these  events  are  within  1000  km. 

GERESS  detections 

The  number  of  detections  (phases)  reported  from  day  092,  1992,  through  day  274,  1992, 
was  43,293,  giving  an  average  of  242  detections  per  processed  day  (178  days  processed). 

Table  3.4.4  shows  daily  and  hourly  distribution  of  detections  for  GERESS. 

Events  automatically  located  by  GERESS 

During  days  092, 1992,  through  274,1992, 3966  local  and  regional  events  were  located  by 
GERESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  22.3  events  per  processed  day  (178  days  processed).  66%  of  these  events  are  within 
300  km,  and  85%  of  these  events  are  within  1000  km. 

J.  Fyen 
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3.5  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  December  1989 
and  was  operated  at  NORSAR  from  1  January  1990  for  automatic  processing  of  data  from 
ARCESS  and  NORESS.  A  new  version  of  IMS  that  accepts  data  from  an  arbitrary  number 
of  arrays  and  single  3-component  stations  was  installed  at  NORSAR  in  October  1991 ,  and 
regular  operation  of  the  system  comprising  analysis  of  data  from  the  4  arrays  ARCESS, 
NORESS,  FINESA  and  GERESS  started  on  15  October  1991.  As  opposed  to  the  first  ver¬ 
sion  of  IMS,  the  one  in  current  operation  also  locates  events  at  teleseismic  distance. 

The  operational  stability  of  IMS  has  been  very  good  during  the  reporting  period.  In  fact 
the  IMS  event  processor  (pipeline)  has  had  no  downtime  of  its  own;  i.e.,  all  data  available 
to  IMS  have  been  processed  by  IMS. 

Events  automatically  located  by  IMS 

During  days  092, 1992,  through  274, 1992, 15,158  events  (local,  regional,  teleseismic) 
were  automatically  located  by  IMS.  This  gives  an  average  of  82.8  events  per  processed 
day  (183  days  processed).  56%  of  these  events  are  within  300  km  of  nearest  station,  and 
69%  of  these  events  are  within  1000  km  of  nearest  station. 

45.8%  of  these  events  were  defined  by  2  regional  phases  and  8.9%  were  defined  by  2  tele¬ 
seismic  phases.  84.8%  of  all  events  had  3  defining  phases  or  less.  18.4%  of  the  available 
detections  (phases)  were  automatically  associated  to  events. 

Events  located  by  analyst  review  of  IMS  results 

During  days  092,  1992,  through  274, 1992, 1 1 ,770  events  (local,  regional  and  teleseismic) 
were  defined  following  analyst  review  of  IMS  results.  This  gives  an  average  of  64  events 
per  processed  day  (184  days  processed).  67%  of  these  events  are  within  300  km  of  nearest 
station,  and  77%  of  these  events  are  within  1000  km  of  nearest  station. 

47.2%  of  these  events  were  defined  by  2  regional  phases  and  4.6%  were  defined  by  2  tele¬ 
seismic  phases.  77.7%  of  all  events  had  3  defining  phases  or  less. 

81.4%  of  the  events  had  regional  phases  only.  15.9%  of  the  events  had  teleseismic  phases 
only. 


Phase  and  event  statistics 

Table  3.5.1  gives  a  sumary  of  phase  detections  and  events  declared  by  IMS.  From  top  to 
bottom  the  table  gives  the  total  number  of  detections  by  the  IMS,  the  number  of  detections 
that  are  associated  with  events  automatically  declared  by  the  IMS,  the  number  of  detec¬ 
tions  that  are  not  associated  with  any  events,  the  number  of  events  automatically  declared 
by  the  IMS,  the  total  number  of  events  defined  by  the  analyst,  and  finally  the  number  of 
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events  accepted  by  the  analyst  without  any  changes  (i.e.,  from  the  set  of  events  automati¬ 
cally  declared  by  the  IMS). 


Apr  92  May  92  Jun  92  Jul  92  Aug  92 

Sep92 

Total 

Phase  detections 

27729 

31278 

34017 

43251 

43331 

37724 

217330 

-  Associated  phases 

8324 

8225 

7702 

8603 

7543 

8673 

49070 

-  Unassociated  phases 

19405 

23053 

26315 

34648 

35788 

29051 

168260 

Events  automatically  declared  by  IMS 

2530 

2493 

2344 

2600 

2346 

2665 

14978 

No.  of  events  defined  by  the  analyst 

2078 

2032 

1802 

2040 

1856 

1854 

11662 

No.  of  events  accepted  without 

1463 

1717 

1580 

1841 

1532 

1503 

9636 

modifications 


Table  3.5.1.  IMS  phase  detections  and  event  summary. 


U.  Baadshaug 
B.  Ferstad 

B.  Paulsen  Gammelby 
B.Kr.  Hokland 
L.B.  Loughran 

3.6  GBF  operation 

Events  automatically  located  by  GBF 

During  days  092, 1992,  through  274,  1992, 11,354  local  and  regional  events  were  located 
by  GBF.  This  gives  an  average  of  62.0  events  per  processed  day  (183  days  processed). 
68%  of  these  events  are  within  300  km  of  nearest  station,  and  89%  of  these  events  are 
within  1000  km  of  the  nearest  station. 

75.9%  of  these  events  were  defined  by  2  regional  phases.  Teleseismic  phases  are  currently 
not  used  by  GBF.  91.1%  of  all  events  had  3  defining  phases  or  less. 

16.4%  of  the  available  detections  (phases,  including  teleseismic)  were  associated  to 
regional  events. 

T.  Kvaerna 
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4  Improvements  and  Modifications 

4.1  NORSAR 

NORSAR  data  acquisition 

No  modification  has  been  made  to  the  NORSAR  data  acquisition  system. 

The  data  are  recorded  on  a  30-hour  circular  disk  buffer  on  the  IBM  system,  and  archived 
onto  1/2  inch  magnetic  tapes.  In  addition  to  this,  the  data  are  now  regularly  transmitted  to 
a  SUN  system  for  recording  on  a  48-hour  circular  disk  buffer. 

NORSAR  detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the  IBM  4381  com¬ 
puter  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2.  In  addition  to  the  detection  processing  done  on 
IBM,  the  dp  program  is  doing  regular  detection  processing  on  a  SUN  system,  using  the 
unix-based  circular  disk  buffer  (see  below).  A  detection  SNR  threshold  of  20.0  triggers 
automatic  saving  of  waveforms  into  CSS  3.0  data  files. 

NORSAR  event  processing 

There  have  been  no  changes  in  the  routine  processing  of  NORSAR  events,  using  the  IBM 
system. 

Routine  event  processing  is  done  on  a  SUN  computer  using  the  “old”  IBM  time  delay  cor¬ 
rection  data  base  that  has  been  converted  to  SUN/UNIX.  The  automatic  solutions  pro¬ 
duced  are  equal  to  or  “better”  than  the  old  system  with  a  lower  false  alarm  rate.  Alert 
messages  are  sent  USGS  for  events  above  magnitude  5.5. 

NORSAR  refurbishment 

As  reported  earlier,  the  main  problem  in  this  refurbishment  is  to  find  24-bit  A/D  convert¬ 
ers  with  small  power  consumption.  Manufacturers  have  so  far  not  been  able  to  present 
production  units  for  evaluation.  One  exception  has  been  Refraction  Technology,  which  has 
allowed  us  to  evaluate  two  prototype  units.  Testing  is  still  going  on,  but  the  units  do  seem 
to  perform  within  specifications.  By  the  writing  of  this  report,  a  new  digitizer  from  Tele¬ 
dyne  Geotech  and  one  from  Science  Horizons  have  been  received  for  evaluation.  All  three 
units  will  now  be  put  out  in  the  field  for  power  evaluation  and  comparison  with  the  NOR- 
ESS  system. 
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4.2  Regional  Arrays 

Detection  processing 

The  routine  detection  processing  of  the  arrays  is  running  satisfactorily  on  each  of  the 
array’s  SUN-3/280  acquisition  systems.  The  same  program  is  used  for  NORSAR,  NOR- 
ESS,  ARCESS,  FINESA,  GERESS,  but  with  different  “recipes”.  The  beam  table  for 
NORESS  and  ARCESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-89/90.  The  beam  table  for 
FINESA  and  GERESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-90/91. 

Detection  statistics  are  given  in  section  3. 

Signal  processing.  Phase  estimation 

This  process  performs  f-k  and  polarization  analysis  for  each  detection  to  determine  phase 
velocity,  azimuth  and  type  of  phase,  and  the  results  are  put  into  the  ORACLE  detection 
data  base  for  use  by  the  IMS. 

Event  Processing.  Plot  and  epicenter  determination 

A  description  of  single-array  event  processing  is  found  in  NORSAR  Sci.  Rep.  No.  2-88/ 
89,  and  NORSAR  Sci.  Rep.  No.  2-89/90. 

New  processes 

The  data  from  the  Apatity  array  are  acquired  on  a  SUN-IPC  system.  Detection  and  signal 
and  event  processing  are  all  performed  on  the  same  computer.  The  new  array  requires 
some  tuning  of  the  recipes  before  data  will  be  routinely  put  into  IMS. 
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5  Maintenance  Activities 

5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  at  the  Maintenance  Center  (NMC)  Hamar,  and 
NDPC  activities  related  to  monitoring  and  control  of  NORSAR,  including  monitoring  of 
NORESS,  ARCESS,  FINESA,  GERESS  and  Poland  (KSP)  partly. 

Activities  involve  preventive  and  corrective  maintenance,  modification  of  equipment, 
installation  planning  (Apatity,  Spitsbeigen),  and  finally  installation  (Apatity). 

NORSAR 

•  Visit  to  subarray  in  connection  with  power  outage 

•  Locating  cables 

•  Replacement  of  SLEM  power  supply 

•  Communication  check  in  connection  with  a  bad  communications  cable 

•  Cable  repair/splicing 

•  Replacement  of  SLEM 

•  Replacement  of  modem 


NMC 

•  Preparing  the  Apatity  and  Spitsbergen  installations 

•  NORSAR  upgrading  activities 

•  Other  tasks  related  to  repair  of  equipment 

NORESS 

•  Cable  location  in  connection  with  road  work 

•  Replacement  of  Hub  card 

ARCESS 

•  Replacement  of  clock,  antenna  and  amplifier 

•  Replacement  of  fiber  optics 

•  Replacement  of  preamplifier  A2  outdoors  and  D2  indoors 

•  Fiber  optic  adjustment 

Apatity 

•  Representatives  from  NTI  (Norwegian  Telecom  International)  and  NORSAR 
installed  NORSAT  B  earth  station  and  SUN  workstation  (May) 

•  Representatives  from  NTI  and  NORSAR  installed  array  equipment  in  the  field  and 
on  the  premises  of  the  Kola  Science  Centre  in  Apatity  (Sept) 
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Subarray/ 

Area 

Task 

Date 

NORSAR 

No  visits  to  the  subarrays 

April 

NMC 

NORSAR  upgrading  and  preparatory  work  related 
to  Apatity  and  Spitsbergen  installations 

April 

ARCESS 

UPS  restart 

April 

NDPC 

Daily  check  of  all  arrays  has  been  done,  i.e.,  NORSAR 
NORESS,  ARCESS,  FINESA,  GERESS  and  Poland 
(KSP, SFP). 

SP/LP  instruments  have  been  calibrated.  FP  and  MP 

April 

were  measured  and  adjusted  (when  feasible  from  NDPC) 
when  outside  specifications. 

NORSAR 

03C 

Cable  pointed  out  for  the  landowner.  Loop  test  in 
connection  with  verification  of  the  communication 

May 

performance. 

02  B 

Visit  in  connection  with  a  power  outage 

25  May 

03C 

Adjustment  of  V,  NS  LP-seismometer 

1  May 

06C 

Adjustment  of  channel  gain  SP  1,2,3  and  5 

Replacement  of  LTA  card  ch  1,2 

18  May 

Apatity 

Representatives  from  NORSAR  and  NTA  installed 
NORSAT  B  earth  station  and  a  SUN  workstation 

May 

NMC 

NORSAR  upgrading 

May 

Preparing  the  Apatity  and  Spitsbergen  installations 

u 

NDPC 

Daily  check  of  all  arrays,  i.e.,  NORSAR,  NORESS, 
ARCESS,  FINESA,  GERESS  and  Poland  (KSP). 
NORSAR  SP/LP  instruments  calibrated  every  week. 

Free  Period  (FP)  and  Mass  Position  (MP)  have  been 
measured  and  adjusted  when  outside  tolerances  (when 
feasible  from  NDPC) 

May 

Table  5.1.  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  including  NDPC 
activities  related  to  the  NORSAR,  NORESS,  ARCESS,  FINESA  and  GERESS  arrays  and 
the  two  3-component  stations  in  Poland  (KSP,  SFP),  1  April  -  30  September  1992. 
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1 

i  Subarray/ 

|  Area 

Task 

Date 

NORSAR 

01B 

Replaced  SLEM  power  supply 

Adjusted  SP  channel  (gain)  2,3,6 

10,11  June 

'  02B 

Replaced  FP  motor 

Adjusted  channel  gain  SP  ch  5,6 

15  June 

01A 

Communication  check  in  connection  with  a  bad  cable 

16  June 

01B 

Cable  SP04  damaged  and  temporarily  repaired 

16  June 

02B 

Adjustment  of  MP/FP  V,  NS,  EW  LP  seismometer 

Also  observed  that  “TEST  ON”  indicator  went  on 
in  connection  with  supposed  spurious  commands 
in  the  ICW  flow  from  NDPC,  cauing  misadjustment 
of  the  6  FP/MP  motors 

19  June 

ARCESS 

1 

Replaced  clock,  antenna  and  amplifier.  Also  replaced 
fiber  optics  B4  A2  outdoors  and  D2  indoors.  Replaced 
DHL70  preamp.  D6  outdoors.  Adjusted  transmit  fiber 
optics  Hub. 

17,18  June 

NMC 

Preparatory  work  in  connection  with  the  Apatity  and 
Spitsbergen  installations 

June 

NDPC 

Daily  check  of  all  arrays,  i.e.,  NORSAR,  NORESS, 

June 

ARCESS,  FINESA,  GERESS  and  Poland  (KSP.SFP  when 
when  running). 

SP/LP  instruments  calibrated.  FP  and  MP  were 
measured  and  adjusted,  and  those  outside  specifications 
were  adjusted  (when  feasible  from  NDPC). 


NORSAR 

02B 

Adjustment  of  channel  gain  SP  1,2,3,4,5,6 
Adjustment  of  channel  gain  LP  V,  NS,  EW 
Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

13  July 

04C 

Adjustment  of  channel  gain  SP  1,2,3,4,5,6 
Adjustment  of  channel  gain  LP  V,  NS,  EW 
Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

14  July 

02B 

Adjustment  of  channel  gain  SP  1,2,3, 4,6 

Adjustment  of  channel  gain  LP  V,  NS,  EW 
Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

17  July 

Table  5.1  (cont) 
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Subarray/ 

Area 

Task 

Date 

02C 

Adjustment  of  channel  gain  SP  1,2,3,4,6 

Adjustment  of  channel  gain  LP  V,  NS,  EW 

Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

20  July 

02B 

Visits  in  connection  with  power  outage  after  lightning 

21,28  July 

03C 

Adjustment  of  channel  gain  SP1,2,3,4,5,6 

Adjustment  of  channel  gain  LP  V,  NS,  EW 

Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

21  July 

NORESS 

Pointed  out  cables  for  the  landowner  in  connection 
with  road  work 

7,8,9,10 

July 

NMC 

Preparatory  work  in  connection  with  the  Apatity  and 
Spitsbergen  installations 

July 

NDPC 

Daily  check  of  all  arrays,  i.e.,  NORSAR,  NORESS,  July 

ARCESS,  FINESA,  GERESS  and  Poland  (KSP). 

Weeks  27, 28  and  29  SP/LP  instruments  have 
been  calibrated.  Free  Period  (FP)  and  Mass  Position 
(MP)  were  measured  and  adjusted  when  outside  tolerances 
(when  possible  from  NORSAR,  Kjeller). 

NORSAR 

02B 

Adjustment  of  MP/FP  V,  NS,  EW  seismometers 

13  Aug 

01B 

Replaced  SLEM  and  modem 

20  Aug 

NORESS 

Replaced  Hub  interface  card 

12  Aug 

NMC 

Continued  installation  planning  in  connection  with 

Apatity  and  Spitsbergen  installations 

Aug 

NDPC 

Daily  check  of  all  arrays  has  been  carried  out,  i.e., 
NORSAR,  NORESS,  ARCESS,  FINESA,  GERESS 
and  Poland  (KSP)  and  Apatity  (partly). 

NORSAR  SP/LP  instruments  calibrated  weeks  27, 28 
and  29.  Free  Period  (FP)  and  Mass  Position  (MP) 
have  been  measured  and  adjusted  when  outside  tolerances 

Aug 

(when  possible  from  NORSAR,  Kjeller). 

Table  5.1  (cont) 
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Subarray/ 

Area 

Task 

bate 

NORSAR 

No  visits  to  the  subarrays 

Sept 

NMC 

Preparations  for  the  installations  in  Apatity  and 

part  of 

Spitsbergen  continued 

Sept 

Apatity 

Representatives  from  NORSAR  and  NTA  satellite 
office  arrived  in  Apatity  22  September,  and  started 
the  installation  the  following  day.  The  work  continued 
through  30  September 

Sept 

NDPC 

Daily  check  of  all  arrays  has  been  carried  out,  i.e., 
NORSAR,  NORESS.ARCESS,  FINESA,  GERESS  and 
Apatity.  NORSAR  SP/LP  instruments  calibrated 
every  week.  Free  Period  (FP)  and  Mass  Position  (MP) 

have  been  measured  and  adjusted  when  outside  tolerances 
(when  possible  from  NORSAR,  Kjeller). 

Table  5.1  (cont.) 


5.2  Array  status 

As  of  30  September  1992  the  following  NORSAR  channels  deviated  from  tolerances: 


01A01 

8  Hz  filter 

02 

8  Hz  filter 

04 

30  dB  attenuation 

02B09 

03C  04 

06C04 

O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Published 

7.1  Extensions  of  the  Northern  Europe  Regional  Array  Network  — 
New  small-aperture  arrays  at  Apatity,  Russia,  and  on  the  Arctic 
island  of  Spitsbergen 

As  of  the  summer  of  1992,  the  Northern  Europe  Regional  Array  Network  contributing 
real-time  data  to  the  NORSAR  Data  Processing  Center  comprised  the  high-frequency 
arrays  NORESS  and  ARCESS  in  Norway,  FINESA  in  Finland  and  GERESS  in  Germany, 
as  well  as  the  two  three-component  stations  KSP  and  SFP  in  Poland.  In  addition,  data 
were  available  in  an  off-line  mode  from  a  three-component  station  installed  in  Apatity  on 
the  Kola  Peninsula  of  Russia  in  June  1991  (see  Mykkeltveit  et  al,  1991).  The  transmission 
of  continuous  data  to  NORSAR  via  satellite  from  the  two  stations  in  Poland  was  discon¬ 
tinued  on  September  30, 1992. 

In  September  and  October  1992,  two  new  small-aperture  arrays  were  added  to  the  net¬ 
work.  These  are  located  near  Apatity,  Russia,  and  on  the  Arctic  island  of  Spitsbergen.  This 
contribution  offers  technical  descriptions  of  these  two  new  arrays  and  their  integration  into 
the  network.  Fig.  7.1.1  shows  the  network  as  of  1  November  1992.  From  approximately  1 
January  1993,  all  data  from  the  six  arrays  of  this  network  will  be  processed  jointly  in  the 
Intelligent  Monitoring  System  (IMS). 

A  new  small-aperture  array  near  Apatity,  Russia 

As  described  in  Mykkeltveit  et  al.  (1991),  a  digital,  high-quality  three-component  station 
was  installed  in  Apatity  on  the  Kola  Peninsula  of  Russia  in  June  1991 ,  as  part  of  an  agree¬ 
ment  of  scientific  cooperation  between  NORSAR  and  the  Kola  Science  Centre  of  the  Rus¬ 
sian  Academy  of  Sciences.  During  June  1991  -  June  1992  data  from  this  station  were 
recorded  locally  in  Apatity  and  copied  to  cartridge  tapes  that  were  sent  to  NORSAR. 

During  the  summer  and  autumn  of  1992,  several  significant  extensions  and  changes  were 
made  to  the  installations  in  Apatity.  In  June,  a  dedicated  satellite  link,  based  on  Norwegian 
Telecom’s  NORSAT  B  system,  was  installed  between  NORSAR  and  the  Kola  Regional 
Seismological  Centre  (KRSC)  in  Apatity.  A  computer-to-computer  Ethernet  connection 
was  established  over  this  link.  This  enables  NORSAR  to  retrieve  continuous  data  from  the 
system  in  Apatity.  At  the  same  time,  several  SUN  workstations  for  data  analysis  were 
installed,  and  a  NOR  AC  array  controller  (Paulsen,  1992)  replaced  the  PC  system  installed 
in  June  1991  with  respect  to  the  data  acquisition  function.  The  configuration  of  the  three- 
component  station  was  otherwise  unchanged.  From  June  1992,  continuous  data  from  the 
Apatity  three-component  station  are  thus  stored  on  Exabyte  cassettes  at  NORSAR. 

In  late  September,  a  small-aperture  array  was  installed  approximately  17  km  to  the  west  of 
KRSC  in  Apatity,  at  the  location  indicated  in  Fig.  7. 1 .2.  The  geometry  of  the  array  is 
shown  in  Fig.  7.1.3.  The  instruments  are  placed  on  two  concentric  rings  plus  one  in  the 
center,  and  the  aperture  is  approximately  1  km.  All  sites  are  equipped  by  a  short  period 
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vertical  seismometer  of  type  Geotech  S-500,  and  the  site  at  the  center  of  the  array  has  in 
addition  two  horizontal  seismometers  of  the  same  type. 

The  field  work  in  preparation  for  the  array  installation  work  was  organized  and  conducted 
by  the  KRSC  and  sponsored  by  the  Kola  Nuclear  Power  Plant  (Kola  NPP).  The  Kola  NPP 
is  located  approximately  30  km  southwest  of  the  array  (see  Fig.  7.1.2).  An  assessment  of 
the  local  seismicity  is  a  topic  of  interest  to  the  Kola  NPP. 

At  each  of  the  seismometer  sites  on  the  A-  and  B -rings,  the  soil  overburden  was  removed, 
and  an  iron  tube  with  a  bottom  plate  was  placed  in  wet  concrete  on  top  of  competent  rock. 
Iron  rods  attached  to  the  underside  of  the  bottom  plate  ensure  the  coupling  of  the  tube  to 
the  concrete  and  underlying  bedrock.  The  vertical  seismometers  were  then  placed  inside 
the  tube,  resting  on  the  bottom  iron  plate.  Power  (commercial  power  is  available  at  the 
central  housing  of  the  array)  is  provided  to  the  seismometers  via  trenched  cables,  and  the 
same  cables  are  used  to  transmit  analog  signals  from  each  seismometer  to  a  centrally 
located  housing  in  the  array.  The  three-component  station  at  the  central  site  AO  is  placed  in 
a  shallow  vault. 

Fig.  7.1.4  shows  in  detail  how  the  data  from  the  new  array  are  acquired  at  the  seismometer 
sites,  digitized  at  the  central  array  site  and  then  transmitted  via  three  radio  channels  to 
Apatity  where  the  NORAC  array  controller  collects  and  timetags  the  data.  The  digitizers 
used  arc  of  types  Nanometrics  three-channel  RD-3  and  six -channel  RD-6,  which  are  both 
16-bit  converters  with  gain  ranging.  Short-period  data  from  the  nine  vertical  sensors  of  the 
array  arc  sampled  at  40  Hz.  The  amplitude  (displacement)  response  function  of  the  short- 
period,  channels  is  given  in  Fig.  7.1.5.  Data  from  the  three  seismometers  at  site  AO  are 
sampled  at  80  Hz,  thus  providing  a  high-frequency  three-component  station  integrated 
with  the  array.  The  seismic  data  from  the  vertical  sensor  at  site  AO  are  thus  sampled  both 
at  40  Hz  (and  used  in  the  processing  together  with  the  vertical  sensors  of  the  A-  and  fi¬ 
rings)  and  at  80  Hz  (and  used  as  part  of  the  high-frequency  three-component  station).  As 
shown  in  Fig.  7.1.4,  there  is  a  certain  redundancy  in  the  data  acquisition,  with  data  from 
sensors  B2  and  B3  being  represented  twice  in  the  data  stream.  This  ensures  that  data  from 
a  useful  array  geometry  are  received  even  in  case  of  failure  of  one  of  the  digitizers  or  one 
of  the  radio  channels.  Also  note  that  timing  signals  in  the  form  of  second  marks  are 
recorded  as  a  separate  time  series,  using  one  of  the  available  channels.  Fig.  7.1.4  also 
shows  the  three-component  broad-band  seismometer  installed  on  the  pier  of  the  basement 
of  the  building  of  the  KRSC  in  Apatity.  The  seismometer  is  of  type  Guralp  CMG-3T,  and 
the  data  are  digitized  using  a  Nanometrics  RD-3  digitizer.  Timing  is  provided  throughout 
the  system  on  the  basis  of  reception  of  GPS  signals,  as  shown  in  the  figure. 

Fig.  7.1.6  shows  the  current  configuration  of  the  data  acquisition  and  analysis  system  at 
the  KRSC.  The  figure  shows  the  local  Ethernet  established,  the  NORAC  array  controller 
that  receives  data  from  the  four  digitizers,  three  SUN  Sparcstations  (kan,  imandra  and 
umb)  and  a  SUN  X-terminal,  and  the  Cisco  router  that  provides  the  gateway  connection  to 
NORSAR  via  the  dedicated  satellite  link.  Also  shown  schematically  is  the  multiplexing 
equipment  used  for  a  phone/fax  connection  via  the  same  satellite  channel.  IDU  is  the  satel¬ 
lite  indoor  unit  containing  the  modem  and  other  communications  equipment.  The  system 
described  in  Fig.  7.1.6  allows  the  staff  at  the  KRSC  to  perform  on-line  processing  as  well 
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as  interactive  analysis  of  the  data  recorded  at  Apatity.  Using  the  satellite  link,  all  data 
recorded  at  NORSAR  from  the  other  arrays  of  the  network  can  also  be  retrieved  by  the 
KRSC  personnel. 

The  overall  picture  related  to  the  Apatity  developments  is  given  in  Fig.  7.1.7.  It  shows  that 
the  data  acquired  in  Apatity  and  retrieved  via  the  NORSAT  B  satellite  link  are  made  avail¬ 
able  to  the  IMS  at  NORSAR. 

A  new  small-aperture  array  on  the  Arctic  island  of  Spitsbergen 

A  small-aperture  array  very  similar  to  the  Apatity  array  described  above  was  installed  on 
the  island  of  Spitsbergen  (see  Fig.  7.1.1)  during  late  October/early  November  1992.  The 
implementation  of  this  array  was  supported  financially  by  Oljeindustriens  Landsforening 
(OLF),  which  is  an  association  of  oil  companies  taking  part  in  oil  exploration  and  produc¬ 
tion  on  the  Norwegian  Continental  Shelf.  The  establishment  and  operation  of  the  commu¬ 
nications  channels  needed  in  order  to  integrate  the  Spitsbergen  array  data  into  IMS  is 
being  supported  by  DARPA  in  the  follow-on  contract  to  the  one  reported  on  here. 

A  suitable  location  for  the  Spitsbergen  array  was  found  during  a  site-selection  survey  in 
August  1991.  The  site  is  located  on  Janssonhaugen  in  Adventdalen  approximately  15  km 
east-southeast  of  Longyearbyen  (see  Fig.  7.1.8),  which  is  the  largest  Norwegian  settlement 
and  also  administrative  center  on  the  island.  Longyearbyen  has  direct  airline  connections 
to  mainland  Norway.  Janssonhaugen  is  a  hill  in  the  middle  of  a  valley  (Adventdalen),  and 
the  array  is  deployed  on  the  plateau  of  this  hill.  The  rocks  at  the  site  are  of  Cretaceous  age, 
covered  by  thin  moraine  of  variable  depth.  The  sensors  are  placed  at  the  bottom  of  6  m 
deep  cased  boreholes.  The  bottom  of  the  boreholes  are  either  in  Cretaceous  rock  or  in 
moraine  material  in  stable  permafrost  conditions  (temperature  approximately  -5°C  all  year 
round  at  a  depth  of  6  m),  such  that  there  is  no  melting/freezing  taking  place  at  this  depth. 

Fig.  7.1.9  shows  the  geometry  of  the  Spitsbergen  array.  All  nine  seismometer  sites  are 
equipped  with  short-period  vertical  seismometers  of  type  Geotech  S-500.  Cables  on  the 
ground  connect  each  of  the  sites  with  a  housing  located  approximately  40  m  south  of  site 
B4.  The  entire  system  with  seismometers,  digitizers,  radio  links  and  array  controller  is 
shown  schematically  in  Fig.  7.1.10.  The  system  design  is  very  similar  to  that  of  the  Apat¬ 
ity  array  described  above.  Power  at  the  Spitsbergen  array  site  is  provided  through  the  use 
of  a  windmill  which  delivers  12V  DC  through  a  battery  bank.  Data  from  the  two  Nano¬ 
metrics  RD-6  digitizers  located  at  the  housing  close  to  site  B4  are  transmitted  over  two 
radio  links  to  Longyearbyen,  where  the  data  are  entered  into  a  NORAC  array  controller 
located  at  a  Norwegian  Telecom  facility.  As  can  be  seen  from  the  figure,  there  is  a  certain 
data  redundancy,  as  data  from  channels  B2  and  B3  are  transmitted  on  both  radio  links.  As 
of  November  1992,  the  seismometers  at  sites  A 1,  A2  and  B3  are  not  yet  connected.  They 
will  be  connected  in  the  spring  of  1993,  as  soon  as  permitted  by  the  prevailing  Arctic  con¬ 
ditions. 

From  Longyearbyen,  the  Spitsbergen  array  data  are  transmitted  via  a  terrestrial  link  to 
Norwegian  Telecom’s  satellite  hub  station  at  Isfjord  Radio,  from  where  a  simplex  64 
kbits/s  satellite  link  (Norwegian  Telecom’s  NORSAT  B  system)  is  used  for  transmission 
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of  the  continuous  data  to  Norway.  The  NORAC  array  controller  in  Longyearbyen  is  also 
connected  to  NORSAR  via  a  dial-up  line,  thus  providing  a  back-link  for  control  and  com¬ 
mand  purposes. 

The  data  from  the  new  Spitsbergen  array  provides  another  data  source  for  the  IMS.  As  this 
array  is  located  in  an  area  of  much  younger  geology  than  the  other  five  arrays,  its  integra¬ 
tion  into  IMS  is  expected  to  present  interesting  seismological  challenges. 


S.  Mykkeltveit  A.  Dahle  J.  Fyen  T.  Kvaerna 

P.W.  Larsen  F .  Paulsen  F.  Ringdal 

I.  Kuzmin,  Kola  Science  Centre 
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Fig.  7.1,1.  The  figure  shows  the  stations  of  the  Northern  Europe  Regional  Array  Network 
as  of  November  1992.  All  stations  of  the  network  are  small-aperture  arrays. 
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Fig.  7.1.2.  Map  of  the  Apatity  region  on  the  Kola  Peninsula  of  Russia,  showing  the  loca¬ 
tion  of  the  Kola  Regional  Seismology  Centre  and  the  location  of  the  new  small-aper¬ 
ture  array  (not  drawn  to  scale;  the  array  aperture  is  1  km).  Also  shown  are  the  mining 
areas  in  Monchegorsk  and  Olenegorsk  to  the  north  of  the  array.  The  seismicity  in  the 
area,  and  in  particular  the  seismicity  associated  with  the  mining  activity  in  the  Khib- 
iny  Mountain  area  to  the  northeast  of  Apatity,  is  described  in  Kremenetskaya  (1991). 
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Fig.  7.0.  Configuration  of  the  new  Apatity  small-aperture  array. 
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Fig.  7.1.5.  Amplitude  (displacement)  response  function  for  the  nine  short-period  channels 
of  the  small-aperture  array  installed  in  Apatity  in  September  1992. 
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Fig.  7.1.7.  Schematic  overview  of  the  field  installations  and  computer  system  in  Apatity  and  the  interfacing  with  the  NORSAR 
Data  Processing  Center  (NDPC)  in  Norway  via  the  NORSAT  B  satellite  link. 
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Fig.  7.1.8.  This  map  of  the  Svalbard  archipelago  with  its  main  island  Spitsbei^en  shows 
the  location  of  the  array  site  at  Janssonhaugen  (J),  the  location  of  the  array  controller 
at  Norwegian  Telecom’s  facility  at  Longyearbyen  (L),  and  the  location  of  the  NOR- 
SAT  B  earth  station  at  Isfjord  Radio  (IR). 


NORSAR  Sci.  Rep.  1-92/93 


November  1992 


Fig.  7.1.9.  Configuration  of  the  new  Spitsbergen  small-aperture  array. 
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chart  for  the  Spitsbergen  array.  1  PPS  (“one  pulse  per  second”)  denotes  timing  signals. 
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7.2  Initial  processing  results  from  the  Apatity  small-aperture  array 

During  the  fall  of  1992,  a  new  regional  array  was  installed  near  the  town  of  Apatity  on  the 
Kola  Peninsula  (see  Section  7.1).  Figs.  7.1.1  and  7.1.2  show  the  location  of  this  array  and 
Fig.  7.1.3  shows  the  array  geometry. 

In  the  following  we  report  on  some  initial  results  from  analyzing  data  from  the  new  array. 
It  is  emphasized  that  these  results  are  preliminary,  since  the  available  data  cover  only  a 
few  weeks.  Thus  a  more  comprehensive  assessment  must  await  the  collection  of  data  over 
a  longer  time  span. 

Noise  spectra 

Fig.  7.2.1a  shows  an  example  of  corrected  noise  spectra  for  the  9  vertical  elements  of  the 
Apatity  array,  taken  at  00.00  hours  (GMT)  on  day  345,  1992.  For  comparison,  the 
ARCESS  spectra  for  the  same  number  of  channels  and  taken  at  the  same  time  are  shown 
in  Fig.  7.2.1b.  From  these  figures,  it  is  seen  that  the  Apatity  array  has  a  higher  noise  level 
than  ARCESS  at  frequencies  above  2  Hz,  whereas  the  noise  levels  are  similar  at  lower  fre¬ 
quencies. 

The  higher  noise  level  at  Apatity  above  2  Hz  is  not  unexpected,  taking  into  account  the 
higher  level  of  industrial  activity  in  this  region.  On  the  other  hand,  the  noise  level  at  the 
Apatity  array  site  is  considerably  lower  than  that  of  the  APA  station  situated  in  the  town  of 
Apatity  (see  NORSAR  Sci.  Rep.  1-91/92). 

Fig.  7.2.2a  shows  24  uncorrected  noise  spectra,  taken  hourly  between  00.00  GMT  and 
23.00  GMT  on  day  344  for  the  Apatity  array  (instrument  AOZ).  For  comparison,  a  similar 
set  of  NORESS  spectra  is  shown  in  Fig.  7.2.2b.  The  range  of  variability  is  similar  for  the 
two  sites,  with  the  highest  noise  levels  corresponding  to  local  daytime. 

Noise  suppression 

Previous  studies  have  shown  that  regional  arrays  are  veiy  effective  in  suppressing  seismic 
noise,  thus  providing  gains  that  are  often  of  the  order  of  Jn,  or  sometimes  even  in  excess 
of  this  value  (N  being  the  number  of  sensors).  Such  better  than  Jn  suppression  occurs 
when  particular  subgeometries  are  chosen,  enhancing  the  suppression  of  noise  at  certain 
frequencies.  As  a  first  check  on  the  capabilities  of  the  Apatity  array  in  this  regard,  we  have 
computed  noise  suppression  curves  and  compared  with  corresponding  results  from  NOR¬ 
ESS. 

We  calculate  an  uncorrected  power  density  spectrum  first  by  prewhitening  60  seconds  of 
data.  Then  we  estimate  the  autocorrelation  for  6  partially  overlapping  windows  (window 
length  12  seconds),  and  compute  power  density  spectra  from  the  average  autocorrelation, 
with  compensation  for  the  prewhitening  filter. 

An  average  spectrum  is  obtained  by  averaging  the  9  individual  channel  spectra  for  the 
array.  The  averaging  is  done  after  a  logarithmic  transform  of  the  spectra,  and  the  standard 
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deviation  for  each  frequency  point  is  calculated.  Each  spectrum  is  pointwise  compared  to 
the  average  spectrum.  If  a  value  is  outside  1.5  standard  deviation  from  the  mean  value,  the 
point  is  considered  an  outlier.  If  a  single  channel  spectrum  has  more  than  60%  outliers,  the 
spectrum  is  excluded  and  a  new  average  spectrum  is  estimated.  A  beam  is  calculated  using 
only  those  channels  for  which  the  spectrum  was  accepted.  The  suppression  is  then  the 
beam  spectrum  divided  by  the  average  spectrum. 

It  should  be  noted  that  calculating  suppression  using  a  beam  with  all  9  elements  compared 
to  using  only  “accepted”  channels  makes  little  difference  in  practice.  When  calculating 
noise  suppression  the  most  stable  results  are  obtained  by  dividing  the  beam  spectrum  with 
the  average  spectrum. 

Fig.  7.2.3a  shows  noise  suppression  in  the  frequency  range  0-20  Hz  for  an  infinite-velocity 
beam  (no  time  delays)  for  the  Apatity  array.  There  are  24  samples,  taken  hourly  on  day 
344, 1992.  For  comparison,  a  corresponding  plot  for  the  9-element  NORESS  B-ring  array 
(AOZ,  A 1-3,  Bl-5)  is  shown  in  Fig.  7.2.3b.  The  Jn  level  is  about  -9.5  dB,  and  the  general 
impression  is  that  the  Apatity  array  is  at  least  as  efficient  in  suppressing  noise  as  the  corre¬ 
sponding  NORESS  subgeometry.  Between  1  and  5  Hz,  the  Apatity  array  has  higher  noise 
suppression  than  the  NORESS  B-ring  geometry,  and  we  attribute  this  to  the  slightly  larger 
Apatity  array  aperture.  Above  5  Hz,  the  two  configurations  are  similar.  At  a  frequency 
near  4  Hz,  there  is  significantly  better  than  Jn  noise  suppression  for  the  Apatity  array.  In 
summary,  it  appears  that  the  spatial  noise  characteristics  at  Apatity  are  similar  to  those 
found  in  other  areas  of  Fennoscandia  and  the  Baltic  Shield. 

Detection  processing 

Since  the  array  was  installed,  the  Apatity  data  have  been  subjected  to  continuous  on-line 
detection  processing  (DP)  at  NORSAR  using  the  detection  algorithm  described  by 
Mykkeltveit  and  Bungum  (1984).  The  initial  beam  deployment  is  shown  in  Table  7.2.1.  At 
an  early  stage,  the  Apatity  detection  lists  were  incorporated  into  the  Generalized  Beam¬ 
forming  (GBF)  process  (Ringdal  and  Kvaema,  1989) ,  which  has  been  in  automatic  opera¬ 
tion  at  NORSAR  for  nearly  3  years.  It  is  noteworthy  that  both  the  DP  and  GBF  could  be 
implemented  for  the  new  array  with  success  even  without  any  special  tailoring  of  parame¬ 
ter  values  and  algorithms. 

As  an  example  of  a  regional  seismic  event  recorded  and  processed  using  Apatity  data,  we 
show  here  results  for  a  presumed  mining  explosion  in  the  Kola  Peninsula  west  of  the  array. 
This  event  is  from  an  area  different  than  the  Khibiny  Massif  where  most  such  explosions 
take  place. 

Fig.  7.2.4  shows  B-ring  seismograms  for  this  event.  Note  the  very  clear  Rg  phase,  which 
indicates  a  shallow  source.  The  detector  output,  and  the  solution  provided  by  the  GBF 
algorithm,  are  shown  in  Table  7.2.2.  A  total  of  10  phases  from  the  four  regional  arrays  in 
Fennoscandia  have  been  associated  with  this  event.  The  estimated  magnitude  (ML)  is 
2.53. 
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Frequency-wavenumber  analysis 

Frequency-wavenumber  solutions  for  the  three  phases  Pg,  Lg  and  Rg  are  shown  in  Fig. 
7.2.5a-c.  In  spite  of  the  small  aperture  of  the  array,  the  peaks  in  the  F-K  diagram  are  well- 
defined:  typically  the  highest  side  lobe  is  at  least  6  dB  below  the  main  peak  in  the  F-K 
plot.  Analysis  of  many  events  has  shown  that  in  particular  the  array  gives  a  very  stable 
determination  of  slowness  and  azimuth  for  the  Rg  phase  in  spite  of  the  broad  peak  in  the 
F-K  diagram.  It  is  remarkable  how  well  an  array  of  only  1  km  diameter  can  resolve  these 
low-frequency  dispersive  waves. 

From  Fig.  7.2.5  we  note  that  the  phase  velocities  of  each  of  the  three  phases  are  consistent 
with  the  phase  type;  the  Rg  phase  having  the  slowest  velocity.  The  estimated  azimuths  are 
also  consistent,  with  a  few  degrees’  difference  only. 

Teleseismic  events 

An  example  of  a  small  teleseismic  recording  is  shown  in  Fig.  7.2.6.  There  is  clearly  a  sig¬ 
nificant  SNR  gain  on  the  array  beam,  and  this  supports  the  previous  statements  on  noise 
suppression.  Naturally,  the  gain  at  teleseismic  signal  frequencies  for  the  Apatity  array  will 
be  lower  than  for  the  full-scale  NORESS,  ARCESS  and  GERESS  regional  arrays.  Never¬ 
theless,  for  areas  of  favorable  signal  focusing  effects,  the  small  Apatity  array  may  still 
show  excellent  teleseismic  detection,  and  many  examples  of  this  have  been  observed. 

In  conclusion,  our  preliminary  analysis  indicates  that  the  Apatity  array  will  be  an  impor¬ 
tant  supplement  to  the  seismic  array  network  in  Northern  Europe,  both  for  regional  and 
teleseismic  event  analysis.  Its  inclusion  into  the  Intelligent  Monitoring  System  (IMS)  will 
in  particular  serve  to  improve  the  location  precision  and  source  characterization  ot  the 
large  number  of  seismic  events  in  the  Kola  Peninsula  and  adjacent  areas. 


F.  Ringdal 
J.  Fyen 
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Table  7.2.1.  Apatity  beam  table,  valid  from  1992-274  (1  October  1992).  The  table  shows 
the  name  of  the  beam,  velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz).  STA/ 
LTA  threshold,  and  configuration.  The  configuration  is  described  with  numDer  of 
sensors  and  a  configuration  code.  Here,  AOAB  means  center  AO  SPZ  plus  A-ring 
plus  B-ring,  and  AOB  means  AO  SPZ  plus  B-ring.  The  designator  (alone)  AO  means 
three-component  80  Hz  data  from  AO.  AH01  -  AH04  are  three-component  horizon¬ 
tal  beams  using  the  AO  high-frequency  system.  A101  -  AI06  are  incoherent  beams 
using  SPZ  channels  only.  (Page  1  of  2) 
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Table  7.2.1  (cont.)  (Page  2  of  2) 
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APA  Single  channel  corrected  spectra. 

60  sec  start  1992-344  23.59.0J.625  Wednesday  09  December 
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Fig,  7.2.1a.  Noise  spectra  corrected  for  system  response  for  the  Apatity  array  for  9  vertical 
channels  at  AO,  the  A-ring  and  the  B-ring.  The  spectra  are  based  on  one  minute  of  data  at 
00.00  hours  GMT  on  day  345, 1992.  The  power  density  is  in  nm2/Hz. 


ARC  Single  channel  corrected  spectra. 

60  sk  start  1992-5*5 00  01  16.556  Itwsdav  10  Ofcembw 


c~-j  t.OE+04. 


a.0E+02. 


S  !.0Et00_ 

* 

o 

Q_ 


1  .OE -02  _ 


0.1  1  10  100 

Frequency  Hz 

Fig.  7.2.1b.  Same  as  Fig.  7.2.1a,  but  for  ARCESS  data  (AO,  A-ring,  B-ring)  taken  at  00.00 
hours  GMT  on  day  345, 1992. 
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APA  hourly  channel  specr  a. 


Sampled  every  how  over  24  hows  Wednesday  09  December 
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Fig.  7.2.2a.  Uncorrected  noise  spectra  for  the  Apatity  array  for  24  one-minute  intervals 
taken  hourly  between  00.00  and  23.00  hours  GMT  on  day  344, 1992.  Each  spectrum 
represents  data  from  the  AO  SPZ  seismometer. 


NRS  hourly  channel  spectra. 


Sampled  every  hour  over  21  hours  dednesdcy  09  December 
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Fig.  7.2.2b.  Same  as  Fig.  7.2.2a  but  for  NORESS  data  (AO  SPZ  seismometer)  taken 
hourly  between  00.00  and  23.00  hours  GMT  on  day  344,  1992. 
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APA  hourly  noise  suppression  spectra. 

Sampled  every  hour  over  24  hours.  Wednesday  09  December 
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Fig.  7.2.3a.  Apatity  array  noise  suppression  by  beamforming  for  the  geometry  comprising 
AO,  A-  and  B-ring  SPZ  sensors.  To  produce  these  curves,  an  infinite-velocity  beam  is 
formed  and  the  spectrum  for  this  beam  is  divided  by  the  average  of  the  single  sensor 
spectra.  The  24  curves  result  from  one  minute  of  data  taken  hourly  between  00.00 
and  23.00  hours  GMT  on  day  344, 1992.  The  horizontal  line  at  -9.5  dB  represents 
Jn  suppression  for  9  sensors. 

NRS  hourly  noise  suppression  spectra. 

bompted  every  how  over  24  hows.  Wednesdcy  09  December 
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Fig.  7.23b.  Same  as  Fig.  7.2.3a  but  for  the  9-element  sub-geometry  of  NORESS  consist¬ 
ing  of  AO,  the  A-  and  B-rings. 
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Fig.  7.2.4.  Plot  of  individual  APA  SPZ  channels  for  the  presumed  mining  explosion  dis¬ 
cussed  in  the  text  (see  Table  7.2.2).  Note  the  very  prominent  Rg  phase,  indicating  a 
shallow  origin. 
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Fig.  7.2.5.  Broadband  F-K  analysis  results  for  the  P,  Lg  and  Rg  phases  of  the  event  shown 
in  Fig.  7.2.4.  The  figure  shows  a)  the  Pg  phase,  b)  the  Lg  phase  and  c)  the  Rg  phase. 
The  three  velocity  circles  correspond  to  8, 6  and  3  km/s.  (Page  1  of  3). 
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1992-298:11.18.41.000  +  1.50  SEC 

WIDE -8AND  SLOWNESS  ANALYSIS 
CONTOURS  IN  DB  BELOW  MAXIMUM  PEAK 


b) 


V,A  3.82  KM/S,  273.68  DEG  65.64  DB 
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Fig.  7.2.5.  (Page  2  of  3) 
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Fig  7.2.5.  (Page  3  of  3) 
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Fig.  7.2.6.  Example  of  teleseismic  event  processing  for  the  Apatity  array.  The  top  six 
traces  are  individual  seismometer  recordings,  filtered  in  the  band  1. 5-3.5  Hz.  The 
two  bottom  traces  are  array  beams,  filtered  and  unfiltered.  Note  the  significant  SNR 
improvement  for  the  filtered  array  beam. 
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7.3  On  the  use  of  regionalized  wave  propagation  characteristics  in 
automatic  global  phase  association 

The  Generalized  BeamForming  (GBF)  method  for  automatic  phase  association  and  event 
location  (Ringdal  and  Kvaema,  1989),  uses  the  philosophy  of  matching  the  predicted 
phase  arrivals  from  hypothetical  events  at  a  predefined  set  of  target  locations  to  the  actual 
detections  at  each  observing  station.  This  approach  gives  us  the  opportunity  of  organizing 
and  using  knowledge  on  regionalized  wave  propagation  characteristics  in  a  simple  and 
well-arranged  way. 

Examples  of  relevant  questions  are: 

•  Which  seismic  phases  are  observable/detectable  at  the  different  stations  from  an 
event  at  a  given  target  location. 

•  What  are  the  expected  travel-times,  slownesses,  and  azimuths  of  these  phases,  and 
the  corresponding  variability. 

•  Wffiat  are  the  amplitude  versus  magnitude  relations  for  the  different  phases,  and  the 
corresponding  variability. 

The  quality  of  the  GBF  results  as  well  as  the  false  alarm  rate  are  dependent  on  how  well 
our  predictions  match  the  actual  observations.  As  illustrated  in  Fig.  7.3.1,  the  detectable 
seismic  phases  may  vary  strongly  from  one  region  to  another,  as  do  the  corresponding 
travel-times.  The  obvious  way  to  obtain  good  predictions  at  a  given  station  is  to  collect 
and  analyze  events  from  a  wide  range  of  epicentral  regions.  For  these  calibration  events, 
independent  bulletin  information  on  hypocenter  location,  origin  time  and  magnitude 
should  be  available. 

Detailed  mapping  of  the  wave  propagation  characteristics  attributed  to  a  given  station  do, 
however,  require  significant  amount  of  data  and  a  substantial  amount  of  work  (Kvaema 
and  Mykkeltveit,  1985;  Mykkeltveit  et  al.,  1990;  Gestermann,  1992;  Sereno  et  al.,  1992). 
The  Intelligent  Monitoring  System  (IMS)  (Bache  et  al.,  1990),  currently  operating  a  net¬ 
work  of  arrays  in  northern  Europe,  has  a  design  that  facilitates  retrieval  of  this  kind  of 
knowledge,  using  a  database  management  system  (DBMS).  During  its  2-3  years  of  opera¬ 
tion  it  has  contributed  significantly  to  improve  the  quality  of  the  automatic  processing 
results  from  this  region  of  the  world  (Bratt  et  al.,  1990). 

The  GSETT-2  database 

Another  interesting  database  is  the  waveforms  and  processing  results  from  the  GSETT-2 
experiment  conducted  in  1991.  During  a  period  of  six  weeks  (April  22  to  June  2),  60  glo¬ 
bally  distributed  stations  (including  single-component  stations,  three-component  stations 
and  arrays)  reported  phase  readings  and  waveforms  to  a  common  database,  see  Fig.  7.3.2. 
At  four  different  data  centers,  automatic  analysis  as  well  as  analyst  review  were  conducted 
to  obtain  precise  event  locations  and  magnitude  estimates. 

In  the  instructions  for  the  conduct  of  the  GSETT-2  experiment  (GSE/CRP/Rev.  4,  1991), 
there  are  detailed  instructions  on  which  parameters  to  use  for  the  preparation  of  the  event 
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bulletin.  This  includes  the  use  of  travel-times  and  slownesses  derived  from  the  Jeffreys- 
Bullen  travel-time  tables.  There  are  also  detailed  instructions  on  which  phases  to  consider, 
as  well  as  their  residual  requirements  on  travel-time  and  slowness  vector  residuals  for  the 
different  types  of  seismic  stations. 

A  stepwise  approach  towards  regionalization 

When  initiating  phase  association  and  event  location  with  the  GBF  method  on  data  from  a 
new  network  like  GSETT-2,  the  natural  starting  point  is  to  use  theoretically  derived  or  glo¬ 
bally  averaged  processing  parameters  like  those  described  in  the  GSETT-2  instructions. 

The  GSETT-2  stations  (see  Fig.  7.3.2)  are  situated  in  very  different  geological  provinces, 
and  may  experience  systematic  deviations  from  the  globally  averaged  travel-times,  slow¬ 
ness  and  azimuth  estimates,  especially  for  events  within  regional  distances. 

A  first  step  towards  regionalization  would  be  to  investigate  if  stations  within  the  same 
geological  province  observe  similar  systematic  anomalies  with  regard  to,  e.g.,  phase 
occurrence,  travel-time,  slowness  or  azimuth,  and  then  accordingly  introduce  corrections 
to  the  parameters  derived  from  the  globally  symmetric  model. 

Next,  this  can  be  refined  by  analyzing  each  station  separately,  and  also,  if  possible,  making 
corrections  for  source  regions  where  the  events  exhibit  anomalous  behavior.  For  example, 
for  seismic  phases  recorded  at  the  GERESS  array  in  Germany,  both  Gestermann  (1992) 
and  Sereno  et  al.  (1992)  have  analyzed  the  characteristics  of  regional  phases  along  differ¬ 
ent  propagation  paths. 

Finally,  for  small  areas  with  recurring  events,  such  as  mines,  nuclear  test  sites  and  seismi- 
cally  very  active  regions,  we  can  use  the  recorded  events  for  calibration  and  obtain  site- 
station  specific  processing  parameters  for  use  by  the  GBF  method. 

In  the  following,  we  will  present  samples  from  results  we  have  obtained  by  searching  the 
GSETT-2  database  that  are  of  relevance  for  improving  the  processing  parameters  for  the 
GBF  method.  In  this  contribution  we  have  used  the  event  bulletin  from  the  Washington 
Experimental  International  Data  Center  (EIDC)  as  the  reference  bulletin,  also  called 
WASCEL  (WAShington  Current  Event  List).  This  bulletin  may  not  contain  totally  inde¬ 
pendent  information  because  the  phases  we  are  to  analyze  have  been  used  to  produce  the 
bulletin.  However,  for  events  with  a  sufficiently  high  number  of  associated  phases,  the 
bias  due  to  this  effect  is  considered  negligible. 

Phase  occurrence 

The  optimum  information  on  the  occurrence  of  seismic  phases  at  a  given  station  would  be 
to  have  a  geographical  map  with  indications  on  which  phases  may  be  observed  from 
events  in  the  different  geographical  areas.  However,  the  derivation  of  such  information 
requires  a  large  number  of  recorded  events,  which  can  only  be  obtained  for  stations  that 
have  been  operational  over  a  long  time  period. 
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If  we  assume  that  the  geology  of  the  region  surrounding  a  given  station  is  rather  uniform, 
the  recording  of  seismic  phases  would  primarily  be  distance-  and  magnitude-  dependent 
(the  effect  of  different  source  functions  is  ignored  at  this  stage).  The  assumption  of  azi¬ 
muthal  symmetry  will  initially  be  made  for  the  derivation  of  wave  propagation  character¬ 
istics  from  the  GSETT-2  database.  We  will  in  the  following  use  phase  labels  reported  by 
the  WASCEL  as  the  reference. 

As  a  first  example  we  show  a  histogram  of  the  occurrence  of  regional  phases  observed  at 
the  ARCESS  array  in  Norway  (ARAO)  as  reported  in  the  WASCEL  (Fig.  7.3.3.a).  It  is 
seen  that  certain  distance  ranges  have  more  reported  phases  than  others,  corresponding  to 
regions  of  high  mining  activity,  like  the  Kola  peninsula  (3-4  degrees),  the  St.  Petersburg 
region  (8-9  degrees),  and  the  Estonian  region  (10-11  degrees). 

This  data  set  indicates  that  beyond  the  critical  distances,  Pn  and  Sn  can  be  expected  to  be 
seen  for  all  regional  distances,  whereas  Pg  is  not  observed  beyond  9  degrees.  Lg  is  only 
observed  within  12  degrees,  and  Rg  only  within  4  degrees. 

It  should  be  emphasized  that  this  data  set  is  sampled  through  a  time  period  of  only  six 
weeks,  and  consequently  only  samples  the  lower  part  of  the  event  magnitude  distribution. 
It  is  therefore  likely  that  phases  from  stronger  events  will  be  observable  at  larger  distances 
than  those  found  in  this  data.  Such  considerations  have  to  be  made  when  constraining  the 
expected  occurrence  of  seismic  phases  at  a  given  station. 

Fig.  7.3.4.a  shows  a  histogram  of  the  occurrence  of  teleseismic  phases  observed  at 
ARCESS.  The  P  and  PKP  phases  generally  follow  the  expected  pattern  of  the  global  travel 
time  curves.  Some  PcP  phases  are  observed,  probably  corresponding  to  events  in  regions 
with  favorable  propagation  paths  and/or  strong  events. 

In  a  similar  way,  such  investigations  can  be  conducted  for  each  of  the  other  GSETT-2  sta¬ 
tions,  but  it  should  be  noted  that  the  GSETT-2  database  probably  contains  too  little  data 
for  a  rigorous  mapping  of  the  occurrence  of  seismic  phases. 

Reported  phase  labels 

Another  interesting  parameter  that  possibly  can  be  utilized  by  the  GBF  algorithm  is  the 
phase  labels  given  by  the  reporting  stations.  For  direct  comparison  with  the  regional  phase 
distribution  of  Fig.  7.3.3.a,  we  show  in  Fig.  7.3.3.b  the  same  type  of  histogram,  but  based 
on  phase  labelling  by  the  Norwegian  National  Data  Center  (NDC).  From  these  figures  it 
looks  like  there  is  good  agreement  with  the  final  phase  labels  given  in  the  WASCEL.  On 
the  other  hand,  from  the  histogram  of  NDC  reported  teleseismic  phases  of  Fig.  7.3.4.b,  we 
find  that  with  very  few  exceptions,  all  phases  were  reported  as  P. 

For  the  most  frequently  observed  phases,  the  relation  between  the  NDC  reported  phase 
labels  and  those  reported  in  the  WASCEL  has  been  studied  in  more  detail.  The  rows  of 
Table  7.3.1  represent  the  ARCESS  phase  labels  reported  by  the  Norwegian  NDC,  and  the 
columns  represent  the  phase  labels  of  the  WASCEL.  This  matrix  table  thus  gives  an  over¬ 
view  of  the  reliability  of  the  NDC  reported  phase  labels. 
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Phase  labels  reported  in  the  WASCEL 


Table  7  3.1. 


Let  us  consider  phases  reported  as  Pn  by  the  Norwegian  NDC: 

186  of  these  phases  were  not  renamed  during  global  phase  association,  4  were  renamed  to 
Pg,4  were  renamed  to  P  and  1  was  renamed  to  be  unknown.  Similar  statistics  are  seen  for 
the  other  phase  types.  In  the  context  of  conducting  event  association  by  the  GBF  method, 
it  is  reasonable  to  introduce  some  general  constraints  on  the  use  of  NDC  reported  phases 
from  the  ARCESS  array.  For  example,  the  following  two  constraints  are  quite  apparent 
from  the  table: 

•  A  phase  reported  to  be  of  P-type  (P,  Pn,  Pg,  PKP,...)  should  not  be  associated  with  a 
hypothetical  arrival  of  S-type.  Out  of  1032  P-type  phases  in  this  data  set,  none  was 
renamed  to  an  S-type  phase. 

•  A  phase  reported  to  be  of  S-type  (S,  Sn,  Lg,  Rg)  should  not  be  associated  with  a 
hypothetical  arrival  of  P-type.  In  this  data  set,  out  of  250  S-type  reports  (Sn,  Lg, 
Rg),  only  one  phase  was  changed  to  a  P-type  phase  by  WASCEL. 

Many  other  constraints  of  more  limited  application  can  also  be  inferred  from  the  table.  For 
example,  a  phase  reported  as  Pn  should  only  be  allowed  to  be  renamed  P  or  Pg  (or  retained 
as  Pn).  On  the  other  hand,  a  phase  reported  as  P  could  be  renamed  to  either  Pn,  Pg,  PP, 
PcP,  PKP,  pP,  sP  or,  of  course,  retained  as  P.  It  is  possible  that  some  more  refinements  may 
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be  made  by  taking  into  consideration  the  difference  in  the  slowness  estimates  between 
core  phases  and  P-phases,  but  this  will  be  the  topic  of  a  separate  investigation. 

Without  going  into  detail,  it  is  clear  that  a  database  of  the  GSETT-2  type  can  provide  a 
number  of  useful  rules  to  constrain  the  reinterpretation  of  phases  reported  by  an  NDC. 
This  in  turn  will  reduce  the  probability  of  false  associations  in  the  global  event  definition 
process. 

Reported  event  locations 

During  the  GSETT-2  experiment,  the  National  Data  Centers  (NDCs)  were  encouraged  to 
provide  information  on  location  parameters  of  events  detected  at  their  participating  sta¬ 
tions. 

These  location  parameters  provide  useful  information  to  the  automatic  phase  association 
procedure.  Depending  on  the  expected  accuracy  of  the  reported  NDC  location,  a  phase 
associated  to  an  event  by  an  NDC,  can  during  the  global  phase  association  process  be 
restricted  to  stem  from  events  in  target  regions  close  to  the  NDC  reported  event  location. 

One-station  event  locations  are  usually  found  by  associating  a  P-  and  an  S -phase  to  a  com¬ 
mon  event,  getting  a  distance  estimate  from  the  S-P  travel-time  difference  and  source 
direction  from  the  azimuth  estimates  of  the  phases.  One-station  event  locations  are  thereby 
most  commonly  obtained  for  regional  events,  where  clear  P-  and  S-phases  are  frequently 
observed.  Reliable  azimuth  estimates  of  both  P-  and  S-phases  are  obtained  by  arrays, 
whereas  three-component  stations  have  the  ability  to  estimate  the  azimuth  of  P-phases. 

To  investigate  the  reliability  of  the  NDC  locations  provided  by  the  regional  arrays, 
ARCESS,  FINESA,  NORESS  and  GERESS,  we  have  in  Fig.  7.3.5  plotted  the  station-to- 
event  distances  of  the  NDC  locations  against  the  station-to-event  distances  of  the  corre¬ 
sponding  locations  of  the  WASCEL.  As  seen  from  the  plots,  there  are,  with  few  excep¬ 
tions,  good  correspondence  between  these  distance  estimates. 

A  few  phases  associated  to  local  or  regional  events  by  the  NDC  have,  however,  been  asso¬ 
ciated  to  teleseismic  events  in  the  WASCEL.  Incidentally,  this  is  a  violation  of  the  rules 
outlined  in  GSE/CRP/Rev.  4  (1991)  (Appendix  B.3.2  (ii):  “Observations  of  phases 
reported  as  originating  from  local  or  regional  distances  may  be  used  only  within  local 
and/or  regional  distances”). 

Conclusions 

We  have  in  this  paper  outlined  a  strategy  for  retrieving  and  organizing  information  on 
regionalized  wave  propagation  characteristics  for  use  by  the  GBF  method  algorithm.  Dif¬ 
ferent  types  of  relevant  information  derived  from  the  GSETT-2  database  have  been  shown 
in  several  examples.  These  examples  illustrated  features  like  the  distance  dependent 
occurrence  of  seismic  phases,  the  reliability  of  reported  phase  labels  and  the  accuracy  of 
one-station  event  locations.  In  future  work  we  plan  to  examine  the  GSETT-2  as  well  as 
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other  data  bases  in  more  detail,  in  order  to  develop  an  optimum  framework  for  application 

of  the  GBF  method  to  a  global  network. 

T.  Kvserna 

U.  Baadshaug 
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B.Okm/$  4.5km/*  3.5km/* 


REDUCED  TRAVEL  TIME  (s) 


Fig.  7  J.1.  Illustration  of  variation  of  relative  importance  of  the  phases  Sn  and  Lg  for  six 
events  with  locations  as  indicated  in  the  map.  The  standard  group  velocities  of  4.5 
and  3.5  km/sec,  commonly  assigned  to  Sn  and  Lg,  respectively,  are  marked  by 
dashed  lines.  The  upper  three  traces  cover  the  distance  interval  from  480  to  550  km, 
while  the  lower  three  traces  correspond  to  epicentral  distances  in  the  range  from 
1,225  to  i,320  km.  The  location  of  the  NORSAR  array  is  denoted  by  a  ring  on  the 
map,  and  the  traces  are  from  NORSAR  seismometer  02B01.  The  data  are  bandpass- 
filtered  1  to  5  Hz.  The  reduction  velocity  is  8.0  km/sec.  From  Mykkeltveii  et  al 
(1990). 
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GSETT-2  stations 


•  Arrays 

A  Three  component 
□  Single  stations 


Fig.  7  J.2.  Map  showing  the  location  of  the  stations  participating  in  the  GSETT-2  experi¬ 
ment. 


NORSAR  Sci.  Rep.  1-92/93 


November  1992 


ARAO  regional  phase  distribution 
4  or  more  defining  phases 


NDC  reported  phase  ids 


Distance  in  degrees 


Fig.  7.3J. 

a)  Histogram  showing  the  distance  dependent  occurrence  of  regional  phases  observed 
at  the  ARCESS  array  in  Norway  (ARAO).  Event  to  station  distances  and  phase 
labels  are  both  taken  from  the  WAShington  Current  Event  List  (WASCEL). 

b)  Histogram  showing  the  ARCESS  regional  phase  labels  reported  by  the  Norwegian 
National  Data  Center  (NDC).  Event  to  station  distances  are  taken  from  the  WASh¬ 
ington  Current  Event  List  (WASCEL). 
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ARAO  teleseismic  phase  distribution 
4  or  more  defining  phases 


Phase  ids  from  WASCEL 


NDC  reported  phase  ids 


Distance  in  degrees 
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Distance  in  degrees 


Fig.  73.4. 

a)  Histogram  showing  the  distance  dependent  occurrence  of  teleseismic  phases 
observed  at  the  ARCESS  array  in  Norway  (ARAO).  Event  to  station  distances  and 
phase  labels  are  both  taken  from  the  WAShington  Current  Event  List  (WASCEL). 

b)  Histogram  showing  the  ARCESS  teleseismic  phase  labels  reported  by  the  Norwe¬ 
gian  National  Data  Center  (NDQ.  Event  to  station  distances  are  taken  from  the 
WAShington  Current  Event  List  (WASCEL). 
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Fig.  73.5.  Plots  showing  the  relation  between  NDC  reported  event  locations  and  those 
given  in  the  WASCEL,  obtained  after  global  phase  association,  for  the  regional 
arrays  ARCESS,  FINESA,  NORESS  and  GERESS.  Two  lines  with  deviations  of +.3 
degrees,  respectively,  are  shown  on  each  plot. 
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7,4  Two  techniques  for  constructing  a  uniform  grid  system  covering 
the  Earth's  surface 

In  many  geophysical  studies  there  is  a  need  to  divide  the  entire  earth's  surface  into  areas  of 
equal  size  and  shape.  An  example,  discussed  in  subsection  7.5  is  the  application  of  the 
Generalized  BeamForming  technique  for  global  seismic  monitoring.  According  to  Rott- 
man  (1960),  there  are  only  five  polyeders  that  can  be  used  for  a  such  division,  and  among 
these,  the  icosaeder  has  the  largest  number  of  subareas  (20  equilateral  triangles).  For  fur¬ 
ther  refinement,  approximations  have  to  be  made. 

We  have  investigated  the  characteristics  of  two  different  techniques  for  quasi-uniform 
gridding  of  the  surface  of  the  globe.  The  first  method  uses  a  reference  meridian  and  a  ref¬ 
erence  latitude  as  a  starting  point,  and  deploys  equidistant  grid  points  along  equidistant 
latitude  circles,  see  Fig.  7.4.1.  To  obtain  complete  and  non-overlapping  coverage,  special 
care  has  to  be  taken  in  the  polar  regions  and  near  the  longitude  opposite  to  the  reference 
meridian. 

The  other  method  is  adopted  from  seismic  prospecting  (Vinje  et  al.,  1992),  where  triangu¬ 
lation  of  the  icosaeder  has  been  used  to  construct  regularly  sampled  wavefronts,  see  Fig. 
7.4.2.  Unlike  the  first  method,  complete  and  non-overlapping  coverage  are  directly 
obtained. 

Properties  of  the  gridding  techniques 

To  visually  compare  the  two  methods,  we  have  in  Fig.  7.4.1  plotted  a  complete  global  grid 
system  produced  by  the  first  technique,  hereafter  called  the  rectangular  method.  In  Fig. 
7.4.2  we  have  plotted  a  similar  grid  produced  by  triangulation  of  the  icosaeder,  hereafter 
called  the  triangular  method.  The  number  of  points  in  these  plots  are  nearly  identical,  i.e., 
645  versus  642.  We  see  that  for  the  rectangular  method  the  shape  of  the  area  spanned  by 
each  grid  point  and  its  eight  neigbours  becomes  heavily  distorted  in  certain  regions.  This 
is  especially  the  case  in  the  polar  regions  and  for  longitudes  far  away  from  the  reference 
meridian.  For  the  triangular  method,  the  shape  of  the  area  spanned  by  each  grid  point  and 
its  six  neighbours  remains  much  more  uniform.  It  should,  however,  be  noticed  that  the  12 
points  constituting  the  original  icosaeder  only  have  five  neighbors. 

To  obtain  complete  coverage  of  the  earth,  we  define  a  circular  region  around  each  grid 
point.  For  the  triangular  grid,  we  have  found  that  circles  with  radii  as  shown  in  Fig.  7.4.4 
completely  cover  the  earth  surface.  We  could  as  well  have  chosen  types  of  geometries 
different  from  the  circle,  like  the  pentagon  or  the  hexagon  where  the  degree  of  overlap  is 
smaller,  but  for  use  in  our  subsequent  analysis,  the  simplicity  of  the  circle  is  preferable. 
For  the  rectangular  grid,  we  have  in  Fig.  7.4.3  plotted  circles  with  the  same  radii  as  in  Fig. 
7.4.4,  and  we  find  that  there  are  certain  areas  that  have  not  been  spanned  by  any  of  the  cir¬ 
cular  regions,  and  hence,  circles  with  larger  radii  are  required.  The  triangular  method 
therefore  gives  more  effective  coverage. 
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Another  appealing  property  of  the  triangular  technique  is  that  each  grid  point  has  a  well- 
defined  pointer  structure  to  its  neighbours.  For  a  given  sub-region,  this  enables  us  to  con¬ 
struct  a  refined  grid  system  in  a  straightforward  manner. 

In  contrast  to  the  rectangular  method,  whore  an  arbitrary  number  of  grid  points  can  be 
specified,  the  number  of  grid  points  obtained  by  the  triangular  method  can  only  take  on 
certain  discrete  values.  These  values  are  given  in  Table  7.4.1,  together  with  the  radii  of  the 
corresponding  circular  regions. 


Number  of  divisions  Number  of  points  Circle  radius  (degrees) 
of  the  icosaeder 


0 

12 

43.8 

1 

42 

21.9 

2 

162 

11.0 

3 

642 

5.5 

4 

2562 

2.7 

5 

10242 

1.4 

6 

40962 

0.7 

Table  7.4.1. 

In  conclusion,  there  are  strong  arguments  for  preferring  the  triangular  method  to  the  rect¬ 
angular  method.  The  only  argument  in  favour  of  the  rectangular  method  is  the  possibility 
of  constructing  a  grid  of  arbitrary  density. 

T.  Kvaerna 
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645  grid  points,  rectangular 


Fig.  7.4.1.  Global  grid  constructed  by  the  rectangular  method.  The  latitude  and  longitude 
steps  were  chosen  to  be  7.95  deg.  The  connection  lines  between  a  grid  point  and  its 
eight  neighbours  are  also  shown. 
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642  grid  points,  triangular 


Fig.  7.4.2.  Global  grid  constructed  by  the  triangular  method.  The  642  points  were 

obtained  by  a  tree-fold  triangulation  of  the  original  icosaeder.  The  connection  lines 
between  a  grid  point  and  its  six  neighbours  are  also  shown. 


100 


NORSAR  Sci.  Rep.  1-92/93 


November  1992 


645  grid  points,  rectangular 


Fig.  7.43.  Circular  regions  of  radius  5.5  deg.  encompassing  each  grid  point  constructed 
by  the  rectangular  method.  Note  that  for  certain  areas  the  coverage  is  not  complete, 
e.g.,  around  70N,  60E  (marked  by  an  arrow). 
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642  grid  points,  triangular 


Fig.  7.4.4.  Circular  regions  of  radius  5.5  deg.  encompassing  each  grid  point  constructed 
by  the  triangular  method.  Note  that  the  coverage  is  complete. 
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7.5  Initial  results  from  global  Generalized  Beamforming 

We  will  in  this  paper  outline  some  of  the  fundamental  concepts  for  using  the  Generalized 
BeamForming  method  (GBF)  to  conduct  phase  association  and  event  location  on  a  global 
scale.  Towards  the  end  of  the  paper  we  will  also  present  initial  results  from  GBF  process¬ 
ing  of  detection  data  from  the  GSETT-2  database.  The  GSETT-2  network  consisted  of 
about  60  globally  distributed  stations  (including  single-component,  three-component  and 
arrays),  see  Fig.  7.3.2.,  which  reported  phase  readings  and  waveforms  to  a  common  data¬ 
base  during  a  period  of  six  weeks.  Detailed  instructions  for  the  conduct  of  the  GSETT-2 
experiment  are  given  in  GSE/CRP190/Rev.4  (1991). 

The  Generalized  Beamforming  method  for  phase  association  and  event  location  (Ringdal 
and  Kvaema,  1989),  has  been  proven  to  work  well  on  a  regional  scale  (distances  within 
approximately  2000  km)  using  continuous  detection  data  from  the  four  regional  arrays 
ARCESS,  FINESA,  GERESS  and  NORESS  (Kvaema,  1990;  Kvaema,  1992).  Taylor  and 
Leonard  (1992)  are  currently  conducting  research  on  applying  this  method  to  teleseismic 
data  from  a  17  station  global  network.  Their  results  indicate  that  the  GBF  method  can  be 
effectively  used  as  a  means  of  conducting  automatic  global  teleseismic  phase  association. 

In  order  to  extend  the  GBF  method  to  work  on  both  regional  and  teleseismic  data  from  a 
large  global  network  we  have  found  it  necessary  to  develop  a  framework  that  facilitates 
testing  of  the  method  with  different  parameter  settings.  We  will  in  the  following  give  a 
detailed  description  of  this  framework. 

Division  of  the  earth  into  regular  target  regions 

The  basic  idea  behind  the  GBF  method  is  to  match  the  predicted  phase  arrivals  from  hypo¬ 
thetical  events  in  a  predefined  set  of  target  regions  to  the  actual  detections  at  each  observ¬ 
ing  station.  We  have  chosen  to  define  these  target  regions  by  their  center  coordinates  and  a 
circular  area  of  a  given  radius.  Initially  we  will  deploy  the  center  of  the  target  regions  at  a 
common  depth  and  ensure  that  the  circular  areas  cover  the  entire  surface  of  the  earth.  In 
practice,  the  target  region  will  constitute  a  volume  around  the  center  grid  point,  and  will 
thus  also  accommodate  hypothetical  events  within  a  given  depth  range. 

In  section  7.4,  we  have  described  two  techniques  for  constructing  a  uniform  grid  system 
covering  the  earth's  surface.  Results  from  this  study  indicate  that  the  method  based  on  tri¬ 
angulation  of  the  icosaeder  (the  triangular  method)  is  preferable.  The  density  of  the  grid 
system  is  given  by  a  single  variable  determining  how  many  times  the  original  icosaeder  is 
divided  by  triangulation.  Figs.  7.5.1  -  7.5.3  show  global  grid  systems  obtained  by  two-, 
three-  and  four-fold  triangulation  of  the  icosaeder.  The  number  of  global  grid  points  in 
these  figures  are  respectively  162, 642  and  2562,  and  the  respective  radii  of  the  corre¬ 
sponding  circular  areas  encompassing  each  grid  point  are  11.0, 5.5  and  2.7  degrees. 
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Parameter  tolerances  due  to  the  areal  extent  of  the  target  region 
Time  tolerances 

On  the  global  map  of  Fig.  7.5.4  we  have  plotted  a  162-point  global  grid  system  with  target 
regions  of  11  degrees  radii.  Let  us  consider  a  P-phase  at  the  NORESS  array  from  a  hypo¬ 
thetical  event  in  the  highlighted  target  region.  As  shown  on  Fig.  7.5.5,  the  P-phase  travel- 
times  from  events  within  this  target  region  can  differ  by  almost  200  seconds.  This  has  to 
be  taken  into  account  in  the  process  of  matching  predicted  and  observed  phase  arrivals, 
and  in  practice,  the  following  procedure  is  used: 

Assume  that  we  have  a  phase  detection  at  the  NORESS  array  with  arrival  time  Tarr  The 
minimum  and  maximum  predicted  travel-times  of  P-phases  from  the  highlighted  target 
region  are  given  by  7Tmjn  and  TT^^.  If  the  observed  arrival  at  NORESS  is  to  be  a  P- 
phase  from  an  event  in  the  target  region,  the  origin  time  of  the  event  would  theoretically 
be  bounded  by  OTmin  and  OTmajc  given  by 


OT  =  T  -  7T 

min  arr  max 

OT  =  T  -  7T 

max  arr  mm 


(1) 

(2) 


Slowness  vector  tolerances 

If  the  observing  station  is  providing  azimuth  and/or  slowness  estimates  of  the  detected 
phases,  these  estimates  can  constrain  the  use  of  the  detections  in  the  phase  association  pro¬ 
cess,  and  thus  prevent  false  associations.  Let  us  again  consider  a  P-phase  at  the  NORESS 
array  from  a  hypothetical  event  in  the  highlighted  target  region  shown  on  Fig.  7.5.4.  To 
cover  the  entire  target  region,  the  P-phases  will  theoretically  span  an  azimuth  and  slow¬ 
ness  range  falling  within  the  solid  curve  of  Fig.  7.5.6.  In  practice,  we  make  an  approxima¬ 
tion  to  the  area  within  the  solid  curve  by  specifying  four  parameters,  i.e.,  the  minimum 
and  maximum  values  of  slowness  and  azimuth.  As  seen  on  Fig.  7.5.6,  this  will  constitute 
an  area  only  marginally  larger  than  the  original.  So  in  order  to  match  the  observed  arrival 
at  NORESS  to  a  P-phase  from  a  hypothetical  event  in  the  target  region,  the  estimated  azi¬ 
muth  and  slowness  would  theoretically  have  to  fall  within  these  bounds. 

Parameter  tolerances  due  to  deviations  from  the  theoretical  model 
Time  tolerances 

In  addition  to  the  parameter  tolerances  compensating  for  the  areal  extent  of  the  source 
region,  we  also  have  to  take  into  account  the  effects  of  sampling  rate,  errors  in  the  estima¬ 
tion  of  arrival-time  and  slowness  vector,  propagation  path,  source  type  and  other  types  of 
random  errors. 

Let  AT  be  the  sampling  rate  in  time  of  the  generalized  beam  for  the  highlighted  target 
region  of  Fig.  7.5.4,  and  let  Tdev-  (early  arrival)  and  Tdev  +  (late  arrival)  be  the  maxi¬ 
mum  allowable  deviations  from  the  predicted  P  travel-time,  taking  into  account  path 
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effects,  onset  time  estimation  errors,  source  effects  and  other  types  of  random  errors.  If  we 
again  consider  the  situation  with  a  phase  detection  at  the  NORESS  array  with  arrival  time 
Tarn  hypothesized  to  be  a  P-phase  from  an  event  in  the  highlighted  target  region,  we  find 
that  the  origin  time  of  the  event  has  to  be  within  the  following  bounds: 

AT 

OTlow  =  Tarr  -  TTmax  -  —  -  Tdev_  (3) 

AT 

OThigh  =  Tarr  -  TTmin  +  -y  +  Tdev+  (4) 

These  are  the  bounds  actually  used  in  the  GBF  matching  procedure.  Initial  residual 
requirements  (Tdev_  and  Tdev+  )  for  P-phases  in  the  GSETT-2  experiment  were  1 .5  sec¬ 
onds,  whereas  for  S-phases  they  were  set  to  7.5  seconds. 

Slowness  vector  tolerances 

One  way  of  including  the  effects  of  random  errors  when  matching  the  estimated  and  pre¬ 
dicted  slowness  vectors,  is  to  require  the  absolute  value  of  the  slowness  vector  residual  to 
be  less  than  a  predefined  value.  In  the  GSETT-2  experiment  different  requirements  were 
used  for  three-component  stations,  high-frequency  arrays  and  short-period  arrays,  reflect¬ 
ing  their  different  abilities  to  correctly  estimate  the  slowness  vector. 

The  procedure  of  comparing  a  slowness  vector  observation  with  the  predicted  slowness 
and  azimuth  range  of  a  target  region,  including  the  residual  requirement,  can  be  described 
by  the  following  example: 

Let  us  again  consider  the  situation  of  a  P-phase  at  the  NORESS  array  from  a  hypothetical 
event  in  the  highlighted  target  region  of  Fig.  7.5.4.  A  phase  detection  at  NORESS  has  a 
slowness  and  azimuth  estimate  given  by  the  asterisk  on  Fig.  7.5.7,  and  according  to  the 
GSETT-2  instructions,  the  corresponding  residual  requirement  is  3.0  sec/deg,  indicated  by 
the  circular  area  encompassing  the  estimated  slowness  vector.  If  there  is  an  overlap 
between  this  circular  area  and  the  sector  segment  corresponding  to  the  expected  range  of 
slowness  vectors  from  the  target  region  (also  see  Fig.  7.5.6),  the  phase  detection  is  said  to 
match  the  P-phase  of  a  hypothetical  event  in  the  actual  taiget  region.  On  the  other  hand,  if 
there  is  no  overlap,  the  observation  is  not  matching  this  hypothesis. 

Additional  constraints  on  the  use  of  a  phase  detection 

When  initiating  the  GBF  procedure  with  a  predefined  set  of  target  regions,  we  compute 
and  store  the  minimum  and  maximum  azimuths  and  distances  to  every  station  in  the  net¬ 
work.  The  availability  of  these  parameters  enables  us  to  constrain  the  use  of  a  phase  detec¬ 
tion  in  a  simple  and  well-organized  way. 

As  an  example,  assume  that  we  have  a  phase  detection  with  peak  frequency  well  above 
10  Hz,  and  that  we  are  trying  to  match  this  detection  to  an  Lg  phase  from  a  hypothetical 
event  in  a  target  region  located  15-20  degrees  away  from  the  recording  station.  Accumu¬ 
lated  statistics  has  shown  that  at  this  station,  Lg  with  peak  frequencies  above  10  Hz  are 
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never  observed  from  events  at  distances  greater  than  5  degrees.  As  such  information  is 
easily  included  in  the  data  structure  of  the  GBF  procedure,  we  can  compare  this  informa¬ 
tion  to  the  actual  distance  to  the  target  region,  and  consequently  find  that  this  phase  detec¬ 
tion  cannot  be  an  Lg  phase  originating  in  this  target  region. 

Another  example:  A  P  and  an  S-phase  recorded  at  an  array  station  have  been  associated 
and  an  initial  event  location  has  been  provided.  Accumulated  event  statistics  at  this  station 
has  given  us  an  estimate  of  the  uncertainty  of  such  one-array  event  locations.  If  a  target 
region  is  located  too  far  away  from  the  one-array  event  location,  these  two  phases  can  be 
excluded  as  originating  from  an  event  in  this  target  region. 

A  third  example:  Several  successive  phase  detections  at  a  given  station  have  been  associ¬ 
ated  with  the  same  event,  either  by  automatic  processing  or  by  an  analyst.  If  we  are  match¬ 
ing  these  detections  to  an  event  in  a  target  region  from  where  Pn  is  expected  to  be  the  first- 
arriving  phase,  only  the  first  of  the  consecutive  detections  can  possibly  be  the  Pn-phase 
(we  ignore  multiple  events  at  this  stage).  Such  context-dependent  information  is  straight¬ 
forward  to  include  in  the  GBF  algorithm. 

These  three  examples  constitute  samples  from  three  different  classes  of  constraints  that 
can  be  imposed  on  the  use  of  a  phase  detection. 

•  Constraints  inferred  from  measurements  on  a  single  phase,  e.g.,  slowness  vector, 
dominant  frequency,  frequency  spectrum,  polarization  attributes,  signal-to-noise 
ratios,  etc. 

•  Constraints  inferred  from  a  reported  event  location  at  a  given  station. 

•  Constraints  inferred  from  the  pattern  of  detections  at  a  given  station.  The  last  two 
types  of  constraints  are  based  on  so-called  context-dependent  information. 

Step-by-step  description  of  the  GBF  method 

Having  defined  some  of  the  key  elements  of  the  GBF  method  for  phase  association  and 
event  location,  we  will  continue  with  a  step-by-step  description  of  the  algorithm. 

•  Definition  of  the  station  network:  Determine  which  stations  to  process  for  phase 
association  and  event  location. 

•  Definition  of  initial  phase  type  candidates:  Decide  the  phase  types  that  may  be 
considered  in  the  phase  matching  process,  e.g.  P,  PKP,  Pn,  Pg,  S,  Sn,  Lg,  Rg,  etc. 

•  Constraining  the  use  of  the  phase  detections:  From  phase  measurements  (domi¬ 
nant  frequency,  slowness,  azimuth,  etc.),  single  station  location  reports  and  other 
single-station  context  information,  we  do,  for  every  phase  detection  at  every  station, 
impose  constraints  on  their  use.  For  each  initial  phase  type  candidate,  every  detec¬ 
tion  is  assigned  a  row  in  a  database  table.  This  row  may  contain  information  on  the 
allowable  distance  range,  depth  range,  and  azimuth  range  of  a  hypothetical  event 
creating  this  hypothetical  phase  type.  If  no  constraints  are  imposed,  the  distance 
range  is  set  to  0-180  degrees,  the  depth  range  to  0-1000  km  and  the  azimuth  range  to 
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+  180  degrees.  If  it  is  determined  that  a  detection  does  not  correspond  to  one  of  the 
initial  phase-type  candidates,  e.g.  Rg,  the  distance  range  of  the  Rg-row  is  set  to  a 
negative  value. 

Construct  a  grid  of  target  regions:  By  using  the  triangular  method,  a  global  grid 
system  of  a  predefined  density  is  constructed,  see  section  7.4. 

Determine  which  of  the  initial  phase  type  candidates  to  consider:  For  every  sta¬ 
tion-target  region  combination,  we  decide  which  of  the  initial  phase  type  candidates 
to  consider  for  phase  matching.  E.g.,  if  a  target  region  is  located  60-70  degrees  away 
from  a  station,  it  makes  no  sense  to  use  regional  phases  like  Pn,  Pg,  Sn,  Lg  or  Rg. 
On  the  other  hand,  if  a  target  region  is  located  at  regional  distances  from  a  station, 
teleseismic  phases  for  that  station  can  be  ignored.  Along  with  the  list  of  phases  to 
consider  from  each  target  region,  we  store  the  expected  travel-times,  slownesses  and 
azimuths,  as  well  as  the  respective  tolerance  limits.  The  decision  of  which  phases  to 
consider  and  their  corresponding  travel-times  and  slowness  vectors  can  be  inferred 
from  general  travel-time  tables  and  great-circle  azimuths,  or,  if  available,  from 
regional  knowledge  on  the  wave  propagation  characteristics  at  a  given  station,  see 
section  7.3. 

Compute  generalized  beams  for  each  target  region:  For  each  target  region,  at  reg¬ 
ular  origin  time  intervals,  we  match  the  observed  phase  detections  in  the  network  to 
the  predicted  phase  arrivals  from  a  hypothetical  event  in  the  region.  The  value  of  the 
generalized  beam  at  a  given  origin  time  is  the  actual  number  of  matching  phase 
detections.  To  avoid  a  list  of  phase  associations  that  may  look  unreasonable  from  a 
seismological  point  of  view,  we  have  imposed  three  constraints  in  the  phase  match¬ 
ing  process.  These  are: 

i)  A  phase  detection  can  be  associated  with  only  one  hypothetical  phase  arrival. 

ii)  A  hypothetical  phase  arrival  can  be  associated  with  only  one  phase  detection, 
preferably  the  one  with  the  smallest  time  residual. 

iii)  The  chronological  order  of  the  associated  phase  detections  have  to  match  the 
chronological  order  of  the  predicted  phase  arrivals.  This  is  to  avoid  inconsis¬ 
tencies  like  an  S-phase  associated  to  arrive  before  a  P. 

Along  with  the  identifications  of  the  matching  detections  (arids),  we  also  store  the 
time  residuals  relative  to  the  center  location  of  the  target  region.  For  stations  provid¬ 
ing  azimuth  and  slowness  estimates,  the  corresponding  residuals  are  also  stored. 

Finding  the  event  location  from  the  best  generalized  beam:  If  our  predictions  of 
phase  arrivals  are  reasonable,  an  event  should  be  reflected  by  a  peak  in  the  general¬ 
ized  beam  representing  the  event  target  region.  However,  like  in  conventional  beam- 
forming,  we  encounter  the  problem  of  sidelobes,  and  will  therefore  have  to  select 
one  generalized  beam  peak  as  representing  the  origin  time  and  location  of  the  event. 
Initially,  we  will  use  the  procedure  of  considering  the  peaks  at  all  generalized  beams 
within  a  predefined  time  interval  and  select  the  one  with  the  highest  number  of 
matching  phases,  and  if  equality,  select  the  one  with  the  smallest  absolute  time  resid¬ 
ual. 


107 


NOR5AR  Sc i.  Rep.  1-92/93 


November  1992 


•  Remove  associated  phases  and  look  for  next  event:  Once  an  event  is  located  and 
the  associated  phases  identified,  these  phases  have  to  be  removed  from  further  con¬ 
sideration.  This  is  done  by  recreating  the  generalized  beams  without  using  the 
phases  associated  to  the  located  event,  and  then  look  for  new  peaks  and  new  events. 

GBF  processing  on  a  global  scale  with  a  large  number  of  target  regions  is  a  computer 
intensive  task.  One  way  of  reducing  the  computational  load  is  to  start  out  with  a  coarse 
grid,  and  then  reprocess  a  denser  grid  surrounding  the  initial  event  location  for  a  time 
interval  around  the  event  origin  time.  In  the  next  sub-section  the  effect  of  the  grid  density 
will  be  illustrated  with  examples  from  GBF  processing  of  GSETT-2  data. 

The  influence  of  grid  density;  an  example  with  an  earthquake  in  Tadjikistan 

As  a  first  test  of  global  Generalized  BeamForming,  we  have  processed  a  1  hour  time  inter¬ 
val  of  GSETT-2  data  that  included  an  earthquake  in  Tadjikistan.  The  event  location  pro¬ 
vided  by  the  Washington  Experimental  International  Data  Center  was: 


Date 

Origin  time 

Latitude 

Longitude 

Depth 

mb 

Ms 

91/05/14 

00:28:45.4 

37.63N 

72.30E 

7  km 

4.3 

4.2 

The  GSETT-2  database  consisted  of  detections  from  60  stations,  and  as  initial  phase  type 
candidates  we  used  the  teleseismic  phases  P,  S  and  PKP,  and  the  regional  phases  Pn,  Pg, 
Sn,  Lg  and  Rg.  This  is  in  accordance  with  the  GSETT-2  instructions  on  the  procedures  for 
preparing  an  initial  event  list.  Jeffereys-Bullen  travel-time  tables  were  used  to  predict  the 
travel-times  and  slownesses,  and  directions  along  the  great-circles  were  used  to  predict 
azimuths.  Regional  phases  were  only  considered  for  distances  within  20  degrees,  and  Rg 
only  within  4  degrees.  The  other  instructions  relevant  to  the  GBF  algorithm,  like  time  and 
slowness  vector  residual  requirements  and  restrictions  on  the  use  of  reported  phases  were 
also  followed  in  the  processing. 

To  investigate  the  influence  of  grid  spacing  on  the  resolution  and  performance  of  the  GBF 
method,  we  first  deployed  a  162 -point  grid  system  covering  the  earth's  surface.  The  corre¬ 
sponding  circular  areas  encompassing  each  grid  point  had  a  radius  of  1 1 .0  degrees.  For  the 
generalized  beams  representing  each  target  region,  we  searched  for  the  maxima  in  a  time 
interval  +  20  minutes  around  the  origin  time  of  the  event.  These  maxima  were  then  inter¬ 
polated  and  contoured  onto  the  global  map  of  Fig.  7.5.8,  where  the  color  scale  represents 
the  number  of  associated  phases  at  the  maxima  of  the  generalized  beams.  A  clear  peak  is 
found  at  target  regions  near  the  event  location,  and  the  overall  maximum  was  13  associ¬ 
ated  phases. 

Fig.  7.5.9  and  7.5.10  represent  the  same  type  of  data,  but  computed  with  denser  grid  spac¬ 
ing,  It  is  clearly  seen  that  the  resolution  is  improved  and  the  sidelobe  effect  is  reduced  with 
denser  grid  spacing.  The  overall  maxima  in  these  figures  were  only  12  associated  phases. 
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indicating  that  with  a  coarse  grid  spacing,  phases  with  large  time  and/or  slowness  vector 
residuals  can  be  associated  as  defining  phases.  However,  by  looking  at  the  results  from 
processing  the  coarse  162  point  grid  system  (Fig.  7.5.8),  we  find  that  it  provides  a  good 
initial  location  of  the  event.  This  location  can  subsequently  be  improved  by  constructing 
and  processing  a  denser  grid  around  the  initial  location. 

We  have  in  this  study  outlined  some  of  the  fundamentals  for  global  Generalized  Beam- 
Forming.  The  procedures  for  processing  target  regions  of  different  areal  extent  has  been 
described,  and  a  step-by-step  description  of  the  procedure  has  been  given.  To  obtain  auto¬ 
matic  event  locations  from  processing  GSETT-2  detection  data,  more  program  coding  is 
needed.  However,  we  have  developed  the  framework  for  GBF  processing  on  a  global  scale 
that  facilitate  research  and  testing  of  the  method.  This  will  become  very  useful  in  the 
development  of  the  complete  system. 


T.  Kvaerna 
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162  grid  points 


Fig.  7.5.1.  A  162-point  global  grid  system  projected  onto  an  azimutal  orthographic  projec¬ 
tion  of  the  earth.  This  grid  system  was  obtained  by  a  two-fold  triangulation  of  the 
icosaeder,  and  each  grid  point  represents  a  target  region  of  1 1  degrees  radius. 
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642  grid  points 


Fig.  7.5.2.  A  642-point  global  grid  system  projected  onto  an  azimutal  orthographic  projec¬ 
tion  of  the  earth.  This  grid  system  was  obtained  by  a  three-fold  triangulation  of  the 
icosaeder,  and  each  grid  point  represents  a  target  region  of  5.5  degrees  radius. 
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2562  grid  points 


Fig.  7 .5. 3.  A  2562-point  global  grid  system  projected  onto  an  azimutal  orthographic  pro¬ 
jection  of  the  earth.  This  grid  system  was  obtained  by  a  four-fold  triangulation  of  the 
icosaeder,  and  each  grid  point  represents  a  target  region  of  2.7  degrees  radius. 
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Distance  37.5  deg  Radius  1 1  deg 


Fig.  7.5.4.  A  162-point  global  grid  system  projected  onto  an  azimutal  orthographic  projec¬ 
tion  of  the  earth.  The  circular  target  regions  are  shown  by  dashed  circles,  and  the  the 
highlighted  target  region  is  given  special  attention  in  this  study.  Its  center  point  is 
located  37.5  degrees  from  the  NORESS  array.  Minimum  and  maximum  azimuthal 
lines  from  NORESS  to  the  target  region  are  also  shown. 
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P-phase  from  37.5  deg. 
Region  radius:  1 1  deg. 
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Fig.  7.5.5.  Figure  showing  the  range  of  predicted  P-wave  travel  times  from  the  highlighted 
target  region  of  Fig.  7.5.4  to  the  NORESS  array. 
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P-phase  from  37.5  deg. 
Region  radius:  1 1  deg. 


Fig.  7.5.6.  Figure  showing  the  range  of  predicted  P-wave  slowness  vectors  at  NORESS 
for  events  in  the  highlighted  target  region  of  Fig.  7.5.4.  The  slowness  vectors  of  the 
circular  boundary  of  the  target  region  are  projeted  onto  the  solid  line  of  this  figure. 
An  approximation  to  the  area  inside  this  solid  line  is  given  by  the  dotted  sector  seg¬ 
ment  defined  by  the  minimum  and  maximum  azimuths  and  the  minimum  and  maxi¬ 
mum  slownesses. 
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Matching  slowness  observation 
Residual  requirement:  3.0  sec/deg. 


Fig.7.5.7.  Figure  illustrating  the  procedure  used  when  matching  observed  and  predicted 
slowness  vectors.  As  explained  in  Fig.  7.5.6,  the  dotted  sector  segment  approximate 
the  expected  range  of  P-wave  slowness  vectors  for  events  in  the  target  region.  The 
slowness  vector  estimate  of  a  detected  phase  is  given  by  the  asterisk,  and  the  maxi¬ 
mum  allowable  slowness  vector  residual  (determined  a  priori)  determines  the  sur¬ 
rounding  circle.  As  seen  on  the  figure,  there  is  an  overlap  between  this  circular  area 
and  the  sector  segment  corresponding  to  the  expected  range  of  slowness  vectors  for 
P-waves  from  the  target  region.  Thereby,  the  phase  detection  is  considered  to  match 
the  slowness  vector  of  the  P-phase  of  a  hypothetical  event  in  the  actual  target  region 
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162  grid  points 
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Fig.  7.5.8.  Contoured  maxima  of  generalized  beams  for  a  162-point  grid  system.  The  max¬ 
ima  of  each  generalized  beam  were  found  by  searching  a  time  interval  of  ±  20  min¬ 
utes  around  the  origin  time  of  an  event  in  Tadjikistan.  The  colour  scale  represent  the 
number  of  associated  phases. 
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642  grid  points 
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Fig.  7.5.9.  Same  as  Fig.  7.5.8,  but  with  a  642-point  grid  system. 
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2562  grid  points 
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Fig.  7.5.10..  Same  as  Fig.  7.5.8,  but  with  a  2562-point  grid  system. 
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7.6  The  Ukrainian  event  of  16  September  1979 

On  28  June  1992,  an  article  appeared  in  the  New  York  Times  of  considerable  interest  to  sci¬ 
entists  working  in  the  field  of  seismic  monitoring: 

Excerpt  from  the  New  York  Times,  28  June  1992 

“(headline)  Izvestiya  Reports  ‘79  A-Test  at  Ukraine  Mine 

MOSCOW,  June  27  (AP)  —  Soviet  scientists  set  off  a  nuclear  blast  in  1979 
next  to  a  Ukrainian  coal  mine,  then  sent  thousands  of  miners  back  to  the 
shaft  a  day  later  without  telling  them,  the  newspaper  Izvestiya  reported. 

The  article,  published  on  Friday,  may  shed  some  light  on  long-standing 
assertions  by  miners  that  a  nuclear  blast  caused  unusually  high  levels  of 
radiation  around  a  town  it  identified  as  Yunokommunarsk. 

Izvestiya  said  officials  have  previously  attributed  the  level  of  radiation  in 
the  area,  which  has  registered  three  or  four  times  normal,  to  industrial 
waste  and  to  the  1986  nuclear  disaster  at  Chernobyl,  625  miles  to  the  north¬ 
west. 

Izvestiya  did  not  report  higher  incidences  of  death,  cancer  or  other  diseases 
in  the  area  near  the  mine,  and  officials  could  not  be  reached  for  comment 
today. 

The  report  said  the  bomb  had  been  detonated  to  see  if  the  explosion  would 
clear  the  mine  of  dangerous  methane  gas.  It  added  that  officials  had  dis¬ 
guised  the  incident  by  staging  a  civil  defense  drill  and  evacuating  the 
town’s  8,000  residents,  most  of  whom  were  miners.” 

It  turned  out  that  the  original  release  from  AP  had  more  information  than  used  by  the  New 
York  Times.  Thus  AP  stated  that  the  detonation  occurred  at  noon  Sept  16, 1979  and  that  the 
yield  was  1/3  kt.  The  location  of  the  town,  Yunokommunarsk,  is  approximately  48.22°N, 
38.30°E. 

Based  on  this  information,  we  examined  available  seismic  event  lists  in  detail  to  see  if  this 
event  could  be  confirmed.  We  searched  for  an  event  with  origin  time  near  09.00.00  GMT, 
which  corresponds  to  noon  Moscow  time.  However,  the  event  was  not  reported,  neither  by 
agencies  using  global  network  data  (ISC,  NEIC)  nor  by  other  available  sources. 

The  Ukrainian  event  was  not  reported  in  the  NORSAR  monthly  bulletin.  However,  the 
NORSAR  automatic  detector  operates  at  a  very  low  threshold,  and  many  small  events, 
especially  at  regional  distances,  are  detected  but  not  included  in  the  monthly  bulletins. 

For  this  reason,  we  decided  to  check  the  original  automatic  detection  lists  in  the  NORSAR 
archive.  It  turned  out  that  the  automatic  (unedited)  NORSAR  bulletin  for  16  Sep  79  con¬ 
tains  an  event  with  the  following  parameters: 
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Origin  time  (GMT):  08.59.53 
Location:  45N,  34E,  Crimea  region 
Magnitude  (mb):  3.3 

The  automatic  plot  for  this  event  has  been  saved  (Fig.  7.6.1),  and  shows  a  high  frequency 
signal,  visible  on  at  least  4  subarrays,  but  with  poor  signal  coherency  across  the  full  NOR¬ 
SAR  array.  The  analyst  did  not  consider  the  event  solution  to  be  of  high  enough  quality  to 
include  it  in  the  final  NORSAR  bulletin,  although  there  is  no  question  that  the  event  is  real 
(and  not  a  noise  detection). 

Unfortunately,  at  the  time  when  this  event  occurred  NORSAR  could  save  only  selected 
data  intervals  on  tape.  At  the  time,  there  seemed  to  be  no  reason  to  save  this  event,  so  the 
digital  data  is  no  longer  available.  Only  the  event  plots  and  the  detector  listings  have  been 
retained. 

The  automatic  location  estimate  (45N  34E)  is  somewhat  different  from  the  location  of 
Yunokommunarsk,  but  well  within  the  uncertainty  for  an  event  with  such  low  signal 
coherency.  The  estimated  origin  time  is  7  seconds  before  the  hour.  As  is  well  known,  the 
traditional  Soviet  PNE  practice  has  been  to  detonate  such  explosions  exactly  on  the  hour. 

Next,  we  tried  to  constrain  the  origin  time  of  the  event,  under  the  assumption  that  its  loca¬ 
tion  was  near  Yunokommunarsk  (see  Fig.  7.6.2).  Under  the  assumption  that  the  event  in 
question  had  an  epicenter  48.22°N,  38.30°E,  its  origin  time  is  estimated  to  be  as  follows: 

08.59.59.6  (using  the  IASPEI 1991  travel-time  tables) 

08.59.56.7  (using  the  Jeffreys-Bullen  tables). 

In  both  cases,  zero  depth  has  been  assumed.  The  signal  arrival  time  was  re-read  visually 
from  the  plot  (the  correction  was  minor  relative  to  the  automatic  arrival  time  estimate). 

As  is  well  known,  the  Jeffreys-Bullen  tables  tend  to  give  2-3  seconds  too  early  origin  time 
estimates  for  Eurasian  explosions,  so  the  data  are  consistent  with  an  origin  time  exactly  on 
the  hour  and  a  location  approximately  as  given  in  the  press  report. 

In  terms  of  nuclear  monitoring  research,  there  are  some  important  lessons  to  be  learned 
from  this  case  study.  Most  importantly,  automatic  detector  information  is  valuable.  It 
should  not  be  discarded,  even  when  the  analyst  determines  that  the  signal  is  of  poor  qual¬ 
ity,  or  the  signal  parameters  are  poorly  defined.  This  is  consistent  with  the  conclusions 
reached  by  the  Geneva  Group  of  Scientific  Experts,  which  has  repeatedly  stressed  that  all 
detected  seismic  signals  should  be  reported  by  the  participating  stations  in  a  future  global 
seismic  network. 

This  case  study  also  highlights  an  important  consideration  for  CTB  verification:  seismic 
methods  will  be  only  one  of  several  monitoring  tools.  If  an  indication  is  obtained  from 
some  other  source  that  a  clandestine  explosion  might  have  taken  place,  it  is  important  to 
be  able  to  confirm  it  through  seismic  recordings.  This  requires  that  at  least  one  station  in  a 
global  network  is  sensitive  enough  to  provide  such  confirmation. 
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Finally,  it  should  be  noted  that  at  the  time  this  event  occurred  there  was  no  advanced 
regional  array  network  in  northern  Europe.  Today,  the  network  of  high-quality  regional 
seismic  arrays  NORESS,  ARCESS,  FINESA,  GERESS,  as  well  as  other  contributing 
arrays,  arc  in  continuous  operation.  Their  recordings  are  regularly  analyzed  using  the 
advanced  Intelligent  Monitoring  System  (IMS).  Had  such  an  event  occurred  today,  it 
would  with  high  probability  have  been  detected,  accurately  located  and  reported.  How¬ 
ever,  an  event  this  small  would  probably  still  be  unidentified  if  the  only  data  available 
were  recorded  at  teleseismic  distances.  Even  with  regional  data  available,  it  can  often  be 
difficult  to  identify  such  low-magnitude  events  reliably  using  current  methods.  Event 
identification  at  low  magnitudes  remains  one  of  the  most  important  research  topics  in  the 
field  of  seismic  monitoring. 

F.  Ringdal,  NORSAR 

P.G.  Richards,  Columbia  University 
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Fig.  7.6.1.  NORSAR  automatic  event  processing  result  for  the  event  described  in  the  text. 
The  upper  three  traces  are  array  beams,  and  the  next  7  traces  are  subarray  beams,  fil¬ 
tered  in  the  band  1. 6-3.6  Hz.  The  vertical  line  marks  the  estimated  arrival  time.  The 
three  lines  of  text  at  the  top  summarize  the  automatically  determined  event  parame¬ 
ters. 
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Fig.  7.6.2.  Map  showing  the  location  of  the  town  Yunokommunarsk  in  Ukraine  (filled 
square)  together  with  the  great  circle  path  to  the  NORSAR  array  (filled  circle).  The 
distance  is  approximately  2200  km. 
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7.7  Induced  seismicity  in  the  Khibiny  Massif  (Kola  Peninsula) 

The  topic  of  this  paper  is  to  review  recent  processes  of  increasing  seismic  activity  in  the 
Khibiny  Massif  in  the  Kola  Peninsula.  It  is  a  typical  example  of  induced  seismicity  caused 
by  rock  deformation  due  to  the  extraction  of  more  than  2-1  tons  of  rock  mass  since  the 
mid-1960s.  The  dependence  of  seismic  activity  on  the  amount  of  extracted  ore  is  demon¬ 
strated.  Some,  but  not  all,  of  the  induced  earthquakes  coincide  with  large  mining  explo¬ 
sions,  thus  indicating  a  trigger  mechanism.  The  largest  earthquake,  which  occurred  on  16 
April  1989  (Ml  =  4.1)  could  be  traced  along  the  surface  for  1200  m  and  observed  to  a 
depth  of  at  least  220  m.  The  maximum  measured  displacement  was  20  cm. 

Introduction 

The  Khibiny  alkaline  Massif  (see  Fig.  7.7.1)  situated  in  the  middle  of  the  Kola  Peninsula 
is  tectonically  unstable.  Its  seismic  activity  has  been  characterized  by  groups  of  a  few 
earthquakes  with  typically  8-10  years  between  groups  (Panasenko,  1969).  The  Massif 
consists  of  different  blocks  separated  from  each  other  by  faults.  The  recent  crustal  move¬ 
ments  are  in  the  range  from  2  to  4  mm/year  (Yakovlev,  1982). 

The  exploitation  of  the  Khibiny  apatite  ores  started  in  1929  and  since  then  about  2.5  •  109 
tons  of  rock  have  been  mined  from  an  area  of  about  10  km2.  This  corresponds  to  a 
decrease  of  the  gravitational  component  by  typically  2.5-3.0  MPa,  and  for  some  parts  of 
the  Massif  by  as  much  as  9-12  MPa. 

At  the  present  time  more  than  108  tons  of  ore  are  extracted  annually  from  three  under¬ 
ground  and  three  open-pit  mines.  The  velocity  of  the  uplift  of  the  near-surface  parts  is  of 
the  order  of  70  mm/year  for  some  tunnels  (Panasenko  and  Yakovlev,  1983). 

Seismic  activity  of  the  Khibiny  Massif 

The  extensive  mining  has  disrupted  the  natural  geodynamic  process  in  the  area,  causing  a 
redistribution  of  crustal  stress,  which  in  turn  has  led  to  increased  seismic  activity.  Fig. 
7.7.2  shows  in  this  respect  the  cumulative  seismic  energy  release  and  the  amount  of 
extracted  ore  for  the  mines  in  the  Khibiny  Massif.  The  seismic  energy  (E)  released  has 
been  calculated  from  the  formula  (Meyer  and  Ahjos,  1985): 

logE  =  12.30+  \.21ML  (1) 

where  E  is  expressed  in  ergs.  The  similarity  of  the  two  curves  strongly  indicates  a  causal 
connection.  It  should  be  noted  here  that  the  earthquake  catalogue  is  considered  homoge¬ 
neous  back  to  1978. 

The  dependence  of  seismic  activity  on  the  extracted  deposit  volume  has  long  been  known 
from  observations  (Glowacka  and  Kijko,  1989).  In  the  Khibiny  Massif  a  more  intensive 
excavation  began  in  the  mid-1960s,  while  the  first  significant  tremors  occurred  in  1981. 
During  that  year  four  felt  earthquakes  (intensity  I  =  3-4  on  the  MSK  scale)  occurred  in  the 
vicinity  of  the  mines.  Some  of  these  earthquakes  were  accompanied  by  sonic  effects. 
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Around  the  same  time,  many  rockbursts  occurred  near  mines  I  and  II  (see  Fig.  7.  /.  1 ),  each 
of  them  displacing  tons  of  rocks  (1-10  m3).  In  fact,  during  five  hours  on  17  May  1981 
more  than  20  rockbursts  occurred.  And  for  the  first  time  the  occurrence  of  an  earthquake  at 
the  time  of  an  explosion  was  observed.  Now  this  situation  is  typical  for  events  in  the  Khib- 
iny  Massif. 

The  largest  earthquakes 

An  earthquake  with  intensity  more  than  5  in  the  nearby  town  of  Kirovsk  occurred  on  29 
August  1982  immediately  after  an  explosion  of  106  tons  in  mine  II  and  at  the  time  of  a 
smaller  quarry  explosion  (4.6  tons)  between  mines  I  and  II. 

The  next  significant  earthquake,  felt  with  intensity  I  =  5  at  Kirovsk,  was  on  19  June  1984. 
Less  than  one  second  before  this  earthquake  there  was  a  small  explosion  (40  kg,  mine  II). 
A  large  block  of  rock  (65x70x70  m3)  was  broken  and  one  part  of  it  was  displaced  relative 
to  another  by  about  5  cm.  After  this  earthquake  there  were  so  many  rockbursts  that  it  was 
impossible  to  carry  out  work  in  this  mine  for  two  weeks. 

The  strongest  such  earthquake  (Ml  =  4.1),  on  16  April  1989,  occurred  almost  simulta¬ 
neously  with  a  large  explosion  (240  tons)  in  mine  I.  This  earthquake  was  felt  with  I  =  8  in 
the  upper  levels  of  the  mine  and  I  =  5-6  at  Kirovsk.  The  maximum  measured  displacement 
was  20  cm,  and  it  occurred  along  a  fault  striking  at  125-135°  and  dipping  at  30-35°  NE. 
The  displacement  was  traced  along  the  surface  for  1200  m  and  observed  to  a  depth  of  at 
least  220  m.  This  event  was  recorded  by  all  seismic  stations  within  a  distance  of  1000  km. 
Fig.  7.7.3  shows  ARCESS  recordings  of  this  earthquake  (distance  420  km). 

The  earthquake  was  accompanied  by  a  great  number  of  aftershocks  attaining  several  hun¬ 
dred  for  the  two  subsequent  months.  On  24  July  1989  another  earthquake  occurred  here. 
Its  size  was  much  less  (see  Table  7.7.1);  however,  remarkable  destruction  again  took  place 
in  this  mine. 

Conclusions 

This  paper  has  discussed  to  principal  topics:  a)  the  question  of  induced  seismicity  caused 
by  rock  excavation  and  b)  the  question  of  mining  explosions  as  a  trigger  for  earthquakes 
and  rockbursts. 

At  present  it  is  too  early  to  make  a  definite  conclusion  about  the  triggering  mechanism 
between  explosions  and  earthquakes.  In  some  cases  pairs  of  such  events  are  separated  by 
some  milliseconds  (for  example,  16  April  1989,  Fig.  7.7.3);  in  other  cases  they  are  sepa¬ 
rated  by  some  seconds  or  more,  sometimes  they  are  not  connected  in  time  (for  example,  2 
June  1992,  Fig.  7.7.4). 
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However,  the  increase  of  seismic  activity  in  the  Khibiny  Massif  caused  by  the  rock  defor¬ 
mation  due  to  the  extraction  of  large  volumes  of  rock  mass  is  clearly  demonstrated. 


E.O.  Kremenetskaya  and  V.M.  Trjapitsin 

Kola  Regional  Seismological  Centre,  Apatity,  Russia 
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No. 

Date 

Origin 

Distance 

Coordinates 

ML 

Seismic 

time 

(km) 

energy (ergs) 

1 

1948 

0923  0000 

17.0 

67.  /ON 

33.60  E 

3.1 

1.58E16 

2 

1955 

0808  172059 

17.0 

67.70  N 

33.60  E 

3.2 

2.24E16 

3 

1955 

0831  2115 

17.0 

67.70  N 

33.60  E 

2.5 

2.80E15 

4 

1960 

0209  210731 

10.0 

67.60  N 

33.60  E 

2.0 

7.08E14 

5 

1974 

0930  091142 

17.0 

67.70  N 

33.70  E 

3.5 

6.31E16 

6 

19/9 

1212  113452 

13.0 

2.1 

9.27E14 

7 

1981 

0416  205634 

200 

2.1 

9.27Ei4 

8 

1981 

3517  075528 

17.0 

2.3 

1.30E15 

9 

1981 

,  .  8  000  /47 

19.2 

2.5 

2.98E15 

10 

1982 

0422  110258 

22.5 

2.1 

9.27E14 

11 

1982 

0829  053335 

17.0 

67.70  N 

33  70  E 

3.3 

3.10E16 

12 

1984 

0619  054731 

17.0 

67.33  N 

33.70  E 

3.9 

1.79E17 

13 

1984 

1030  105158 

35.0 

2.2 

1.24E15 

14 

1984 

1030  142148 

17.5 

67.68  N 

33.72  E 

2.3 

1.66E15 

15 

1987 

0725  161339 

22.5 

67.66  N 

33.90  E 

2.9 

9.60E15 

16 

1988 

0113025153 

25.0 

67.73  N' 

33.83  E 

2.6 

4.00E15 

P 

1988 

0118  020948 

24.0 

67  65  N 

33.96  E 

2.6 

4.00E15 

18 

1988 

0120  121510 

6.0 

67.60  N 

33.51.  E 

26 

4.00E15 

19 

1988 

0211  124113 

27.0 

2.2 

1.20E15 

20 

1988 

0304  231701 

18.0 

67.70  N 

33.70  E 

2.1 

9.30E14 

21 

1988 

0416  115725 

18.0 

67.66  N 

33.75  K 

2.1 

9.30E14 

22 

1988 

0622  013408 

10.5 

67.65  N 

33.47  E 

2.4 

2.20E15 

23 

1988 

1006  094741 

32.0 

67.61  N 

34.19  E 

3.3 

3.1  OH  1 6 

24 

1988 

1123  211108 

23.0 

67.60  N 

33.80  E 

2.5 

3.00E15 

25 

1989 

0203  102741 

33.7 

67.80  N 

33.90  E 

2.2 

1.20E15 

26 

1989 

0416  063442 

18.0 

67.61  N 

33.81  E 

4.1 

3.20E17 

27 

1989 

0707  114924 

26.5 

67.71  N 

33.93  E 

3.4 

4.10E16 

28 

1989 

0724  223234 

15.0 

67.60  N 

33.78  E 

2.5 

3.00E15 

29 

1989 

0804  012618 

21.0 

67.60  N 

33.90  E 

2.1 

9.30E14 

30 

1990 

0210  163907 

37.1 

67.89  N 

33.48  E 

2.2 

1.20E15 

31 

1990 

0403  075417 

26.5 

67.60  N 

33.90  E 

2.1 

9.30E14 

32 

1990 

0612  113532 

36.7 

67.63  N 

34.29  E 

2.1 

9.30E14 

33 

1990 

0621  130954 

20.5 

67.67  N 

33.83  E 

2.3 

1.70E15 

34 

1990 

0621  231904 

15.5 

67.67  N 

33.66  E 

2.2 

1.20E15 

35 

1990 

0622  015115 

16.5 

67.61  N 

33.81  E 

2.4 

2.20E15 

36 

1990 

0624  065848 

19.4 

67.61  N 

33.88  E 

2.4 

2  20E15 

37 

1990 

0625  062215 

16.2 

67.56  N 

33.82  E 

2.3 

1.70E15 

38 

1990 

0625  113116 

16.1 

67.61  N 

33.80  E 

2.2 

1.20E15 

39 

1990 

0625  215128 

15.3 

67.67  N 

33.66  E 

2.3 

1.70E15 

40 

1990 

0626  052403 

19.1 

67.72  N 

33.30  E 

2.4 

2.20E15 

41 

1990 

0627  004129 

24.4 

67.77  N 

33.49  E 

2.3 

1  JOE  15 

42 

1990 

0627  025428 

19.6 

67.67  N 

33.79  E 

2.2 

1.20E15 

43 

1990 

0627  051520 

23.6 

67.59  N 

33.99  E 

2.4 

2.20E15 

44 

1990 

0627  064141 

29.0 

67.66  N 

34.07  E 

2.8 

7.20E15 

45 

1990 

0628  015919 

22.8 

67.74  N 

33.22  E 

2.6 

4.00E15 

46 

1990 

0629  041211 

19.4 

67.69  N 

33.75  E 

2.8 

7.20E15 

47 

1990 

0630  045139 

20.4 

67.62  N 

33.90  E 

2.4 

2.20E15 

48 

1990 

0630  064224 

28.0 

67.64  N 

34.10  E 

2.6 

4.00E15 

49 

1991 

0305  000836 

23.7 

67.55  N 

34.00  E 

2.1 

9.30E14 

50 

1991 

0505  054802 

28.5 

67.65  N 

34.07  E 

2.1 

9.30E14 

51 

1991 

0815  083545 

24.1 

67.63  N 

33.98  E 

2.1 

9.30E14 

52 

1991 

1229  151425 

27.0 

67.70  N 

33.96  E 

2.1 

9.30E14 

53 

1991 

1229  152044 

25.0 

67.57  N 

34.03  E 

2.1 

9.30E14 

Table  7.7.1.  Known  earthquakes  (ML>2.0)  in  the  Khibiry  Massif  1948-1991.  The  dis¬ 
tance  to  the  APA  station  in  Apatity  is  indicated. 
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33.00  33.25  33.50  33.75  34.00  34.25  34-50 


Fig.  7.7.1.  The  position  of  mines  (I- VI)  in  the  Khibiny  Massif  together  with  fault  struc¬ 
tures  and  earthquakes.  Mines  I,  II  and  III  are  underground,  whereas  IV,  V  and  VI  are 
open-pit  mines.  The  location  of  the  Apatity  seismic  station  (APA)  is  also  shown. 
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Fig.  7.7.2.  Cumulative  distributions  versus  time  of  extracted  masses  of  the  apatite  mineral 
(1)  and  BeniofT’s  graph  of  released  energy  (2)  for  the  Khibiny  Massif.  Note  the  sim 
ilarity  between  the  two  graphs. 


NORSAR  Sci.  Rep.  1-92/93 


November  1992 


NORSARSci.Rep.  1-92/93 


November  1992 


